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ABSTRACT 


An  analytical  model  has  been  developed  which  will  predict  the 
performance  and  pollutant  emissions  of  gas  turbine  combustors.  The 
entire  gas  turbine  combustor  Is  approximated  as  a  collection  of  per¬ 
fectly  stirred  zones.  Within  each  zone  a  general  hydrocarbon  contus- 
tlon  mechanism  Is  used  to  predict  the  gas  composition  and  temperature. 
The  zone  volumes  and  sizes  are  assigned  from  consideration  of  the 
theoretically  predicted  gas  flows  thereby  approximating  the  mixing 
behavior  of  the  system.  The  resulting  combustor  model  Is  seen  to 
possess  several  Intrinsic  merits,  namely: 

1.  It  is  simple  enough  to  facilitate  solutions; 

2.  both  mixing  and  kinetics  are  considered; 

3.  there  Is  an  exact  correspondence  between  the 
flow  zones  In  the  model  and  In  a  combustor. 

Selected  predictions  of  the  overall  model  ror  a  "typical" 
aircraft  combustor  are  presented.  These  results  are  seen  to  be 
qualitatively  accurate  and  fall  1  r.  the  range  of  values  typically  ob¬ 
served  In  practical  systems. 
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CHAPTER  I 


INTRODUCTION  AND  OBJECTIVES 

It  Is  now  apparent  that  the  gas  turbine  engine  will  take  an 
Important  part  in  the  field  of  vehicular  propulsion  in  the  near  future. 
These  engines  are  being  considered  as  alternative  powerplants  both 
for  heavy  vehicles  (trucks  and  military  vehicles}  and  for  automobiles. 
For  vehicular  applications  a  regenerator  is  usually  included  in  the 
engine  cycle,  and  this  regenerator  causes  an  Increase  in  the  combustor 
Inlet  air  temperature.  A  similar  but  less  severe  inlet  temperature 
Increase  is  encountered  in  aircraft  engines  as  the  engine  pressure 
ratios  and  flight  Mach  numbers  increase.  Lifting  engines  for  STOL/VTOL 
applications  also  experience  periods  of  hot  gas  ingestion  which  raises 
the  inlet  temperature  significantly.  For  all  of  these  situations  large 
quantities  of  empirical  design  data  are  not  available  at  elevated 
combustor  inlet  temperatures.  Rather  than  attempting  to  scale  up 
low  Inlet  temperature  data,  it  is  felt  that  an  analytical  model  of 
the  combustor,  valid  for  all  inlet  temperatures,  will  be  eminently 
more  useful. 

The  current  enphasis  on  reducing  the  emission  of  atmospheric 
pollutants  from  all  types  of  engines  has  3lso  focused  attention  on 
the  utility  of  the  analytical  modelling  of  combustion  systems.  In¬ 
creased  application  of  gas  turbine  engines  can  only  Increase  interest 
in  the  control  of  their  emissions.  The  primary  gas-phase  pollutants 
emitted  from  gas  turbines  are  carbon  monoxide  and  nitric  oxide.  The 
emissions  of  unburned  hydrocarbons  are  usually  negligible,  but  smoke 
emissions  can  be  significant.  The  '-evels  of  carbon  monoxide  and 
nitric  oxide  are  clearly  dependent  on  the  kinetic  and  mixing  pro¬ 
cesses  occuring  Inside  the  combustor.  Smoke  emissions  are  controlled 
primarily  by  the  concentrations,  temperatures,  and  flow  pattern  in 
the  primary  combustion  zone.  All  of  the  processes  mentioned  can  be 
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incorporated  into  an  analytical  model  to  permit  prediction  of  emission 
levels. 

The  development  of  the  overall  gas  turbine  combustor  model 
followed  the  same  logical  sequence  as  will  be  followed  here.  Firstly, 
the  kinetics  of  hydrocarbon  combustion  are  investigated  results. g  in 
a  general  hydrocarbon  combustion  theory  for  a  single  perfectly  stirred 
reactor.  The  kinetics  of  nitric  oxide  forratlon  are  incorporated  into 
the  combustion  mechanism.  Secondly,  several  models  of  the  turbulent 
mixing  phenomena  are  considered.  Thirdly,  the  flow  processes  Inside 
the  combustor  are  examined.  Combining  these  three  studies,  a  complete 
analytical  model  is  postulated  and  used  to  make  quantitative  predi- 
tions. 
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CHAPTER  II 

A  GENERALIZED  HYDROCARBON  COMBUSTION 
ANALYSIS  FOR  A  PERFECTLY  STIRRED  REACTOR 

The  overall  gas  turbine  combustor  model  to  be  used  in  this  work 
was  presented  by  Hammond  and  Mellor  (1970a)  and  Hammond  and  Mellor 
(1970b).  Certain  modifications  have  been  made  to  the  original  model, 
but  it  still  remains  basically  a  combination  of  perfectly  stirred  re¬ 
actors.  A  kinetic  analysis  for  the  combustion  of  propane  in  a  per¬ 
fectly  stirred  reactor  was  previously  developed  (Hammond  and  Mellor, 
1970a, b)  for  use  in  this  model.  This  original  analysis,  upon  detailed 
Investigation,  was  shown  to  possess  severa1  undesirable  limitations 
which  subsequently  resulted  in  its  discard.  Firstly,  the  analysis 
was  valid  only  for  propane  which  is  an  unrealistic  fuel  from  a  pract¬ 
ical  viewpoint.  Secondly,  an  ill  defined  empirical  parameter  was 
Involved  in  the  analysis  (cf.  Schneider,  1960).  Thirdly,  predictions 
could  only  be  obtained  for  a  very  limited  range  of  equivalence  ratios 
(0.8  ±<t>  <_  1.2).  Finally,  the  limitation  on  equivalence  ratio  pre¬ 
cluded  any  realistic  evaluation  of  accuracy  by  comparison  with  exper¬ 
imental  data. 

The  analysis  which  follows  was  developed  to  allow  a  priori  pre¬ 
dictions  of  the  combustion  performance  for  a  general  hydrocarbon/air 
system  in  a  single  stage  perfectly  stirred  reactor.  All  of  the  cited 
limitations  inherent  in  the  previous  analysis  have  been  successfully 
removed . 

A.  The  generalized  hydrocarbon  combustion  mechanism 

The  exact  kinetic  mechanisms  for  hydrocarbon/air  combustion  are 
poorly  understood  except  for  the  simplest  hydrocarbon  species.  In 
those  cases  where  "complete"  mechanisms  have  been  postulated,  the 
mechanisms  contain  a  nunber  of  elementary  reactions  for  which  accurate 


rate  data  are  essentially  unavailable.  There  have  been  many  attempts 
to  model  the  combustion  process  in  a  global  manner;  a  number  of  these 
are  summarized  in  Hammond  and  Mellor  (1970a)  and  in  Edelman  and  Fortune 
(1969).  Our  intent  here  is  to  formulate  an  approximate  mechanism 
which  will  apply  to  a  general  hydrocarbon,  and  in  order  to  do  this  we 
will  first  examine  some  of  the  "complete"  mechanisms. 

Chinitz  (1965)  used  the  methane  combustion  mechanism  shown  in 
Table  1  to  study  the  constant-pressure  combustion  of  methane/air  systems. 
His  predicted  temperature-time  histories  compared  very  well  with  ex¬ 
perimental  values,  but  his  predicted  species  concentration-time  his¬ 
tories  deviated  significantly  from  experimental  values.  The  predic¬ 
tions  of  ignition  delay  times  were  about  a  factor  of  100  lower  than 
experiment;  however,  Chinitz  felt  that  mixing  limitations  in  the 
physical  experiment  can  be  charged  with  some  of  the  discrepancies. 

Note  that  Reactions  [1-1]  through  [1-5]  constitute  a  process  of  methane 
partial  oxidation  to  carbon  monoxide  and  water;  that  Reaction  [1-6] 
is  the  only  carbon  monoxide  destruction  reaction;  and  that  the  remainder 
of  the  mechanism  is  concerned  only  with  free  radical  production. 

Morgan  (1965)  proposed  the  mechanism  shown  in  Table  2  for  methane/ 
air  combustion  in  a  perfectly  stirred  reactor.  Using  this  mechanism 
he  was  able  to  derive  certain  global  rate  expressions  which  agreed 
reasonably  well  with  his  own  data.  This  mechanism  differs  only 
slightly  from  that  of  Chinitz,  and  can  be  similarily  interpreted  as 
consisting  of  partial  oxidation  to  CO,  carbon  monoxide  destruction, 
and  free  radical  production  processes.  Morgan  views  Reactions  [2-1] 
and  [2-10]  as  being  rate-limiting. 

The  mechanism  shown  in  Table  3  was  developed  by  Bowman  and  Seery 
(1968)  and  also  by  Bowman  (1968)  from  shock  tube  studies  of  methane 
oxidation.  Both  works  showed  reasonable  agreement  between  predicted 
and  experimental  values  of  the  concentration-time  histories  and  ignition 
delay  time  (defined  as  the  time  to  the  rapid  increase  in  pressure)  for 
mixtures  around  stoichiometric  and  temperatures  above  2000°K.  Bowman 
and  Seery  attributed  the  failure  to  predict  behavior  over  the  full 
range  of  terrperature  and  stoic  hi  one  try  to  poorly  understood  reactions 
involving  HCO  which  were  omitted  from  their  .mechanism. 
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Table  1 . 

Methane  combustion  mechanism  of  Chinltz  (1965) 


Reaction  Ho. 

ch4  tOH  :  ch3  +  H2Q  [1-1] 

CH4  +  0  t  CH3  +  OH  [1-2] 

CH3  +  02  t  HCHO  +  OH  [1-3] 

HCHO  +  OH  J  HCO  +  H20  [1-4] 

HCO  +  OH  t  CO  +  H20  [1-5] 

CO  +  OH  t  CO,  +  H  [1-6] 

H  +  H20  t  OH  +  H2  [1-7] 

0  +  H2  ♦  OH  +  H  [1-8] 

OH  +  OH  *  H2°  +  O  [1-9] 

H  +  0,  j  OH  +  0  [1-10] 

0+0+M  J02+M  [1-11] 

t  H2  +  M  [1-12] 


H  +  H  +  M 
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Table  2. 

ffethane  combustion  mechanism  of  Morgan  (1967) 


Reaction 

CH4  +  OH  t  ch3  +  H20  C2"1! 

ch4  +  H  t  ch3  +  h2  C2'23 

CH3  <  02  t  HCHO  OH  [2-33 

ch3  +  0  t  ch2o  *  H  C2"43 

HCHO  +  OH  t  HCO  +  H20  C2"53 

HCHO  +  M  t  HCO  +  H2  C2'63 

HCO  +  OH  t  CO  ♦  h2o  C2"73 

HCO  +  CH3  *  co  +  ch4  C2_83 

HCO  +  H  t  CO  +  H2  [2-93 

CO  +  OH  i  C02  +  H  [2-10] 

H  +  02  t  OH  +  0  [2“1 1 3 

H20  +  0  t  OH  +  OH  [2-12] 

H  +  H20  t  OH  +  H2  [2-13] 
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Table  3. 

Methane  combustion  mechanism  of  Bowman  and  Seery  (1968) 


Reacti on 

Ho. 

CH4  +  M 

-4 

<4* 

CH3  +  H  +  M 

[3-1] 

ch4  +  o2 

t 

ch3  +  ko2 

[3-2] 

02  +  M 

t 

0  +  0  +  M 

[3-3] 

ch4  +  0 

t 

CH3  +  OH 

[3-4] 

CH4  +  H 

t 

ch3  +  h2 

[3-5] 

CH4  +  OH 

* 

CH3  +  H20 

[3-6] 

ch3  +  0 

-► 

4- 

HCHO  +  H 

[3-7] 

CH3  +  02 

* 

HCHO  +  OH 

[3-8] 

HCHO  +  CH 

-♦ 

4- 

MCO  +  H20 

[3-9] 

HCO  +  OH 

-► 

4- 

CO  +  h2o 

[3-10] 

CO  +  OH 

t 

C02  +  H 

[3-11] 

H  +  02 

-4 

4- 

0  +  OH 

[3-12] 

o  +  h2 

t 

H  +  OH 

[3-13] 

o  +  h2o 

“► 

4- 

OH  +  OH 

[3-14] 

h  +  h2o 

* 

H2  +  OH 

[3-15] 

H  +  OH  +  M 

* 

H20  +  H 

[3-16] 

CH3  +  02 

-4 

4- 

HCO  +  H?0 

[3-17] 

HCO  ♦  02 

-+■ 

4- 

CO  ♦  ho2 

[3-18] 

H  +  H02 

t 

OH  +  OH 

[3-19] 

HCO  +  M 

-4- 

4- 

H  ♦  CO  +  M 

[3-20] 
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Although  the  methane  combustion  mechanism  of  Bownan  and  Seery 
Is  rather  more  complex  than  the  previous  two  mechanisms.  It  still 
possesses  the  three  rather  well  defined  processes.  These  processes 
are  partial  oxidation  (Reactions  [3-1],  [3-2],  [3-4]-[3-10],  [3-17]- 
[3-20]),  carbon  monoxide  destruction,&nd  free  radical  production.  The 
parametric  studies  of  Bowman  and  Seery  (1968)  indicated  the  dominant 
reactions  to  be  [3-11],  [3-4]-[3-6],  and  [3-10]-[3-17], 

Bowman  and  Seery  (1968)  also  developed  two  mechanisms  for  acety¬ 
lene  combustion  (cf.  Table  4);  the  two  mechanisms  differ  only  in  the 
choice  of  the  oxygen  atom  attack  reaction  on  acetylene  from  Reaction 
[4-2a]  or  [4-2b].  The  ignition  delays  predicted  using  each  mechanism 
were  In  good  agreement  with  experimental  values  and  with  each  other. 

The  predicted  temperature-time  and  concentration-time  histories  were 
different  for  the  two  mechanisms,  and  the  predicted  discrepancies 
were  dependent  on  the  rates  of  the  reactions  involving  acetylene  and 
acetylene  fragments. 

Surprisingly,  this  acetylene  combustion  mechanism  can  also  be 
divided  into  the  processes  of  partial  oxidation  (Reactions  [4-ll]-[4-8]), 
carbon  monoxide  destruction  (Reactions  [4-9]),  and  free  radical  pro* 
duction.  Since  acelylene  has  a  triple  bond  structure  and  methane  only 
single  bonds,  one  might  expect  a  significantly  different  combustion 
mechani sm. 

A  detailed  discussion  of  possible  propane  combustion  mechanisms 
was  presented  by  Hammond  and  Mel lor  (1970a),  All  of  the  mechanisms 
considered  In  that  work  could  also  be  divided  Into  the  aforementioned 
three  processes.  Further,  the  partial  oxidation  process  was  assumed 
to  be  Infinitely  fast  compared  to  the  rate  limiting  process  of  carbon 
monoxide  destruction  through: 

CO  +  OH  -  C02  +  H  [24] 

which  was  the  only  destruction  reaction  (in  agreement  with  all  of  the 
above  mechanisms).  The  two  body  branching  reactions: 

OH  +  H2  t  H20  +  H 

H  +  o2  i 


OH  +  0 


[2] 

[3] 


Table  4 


The  acetylene  combustion  mechanisms  of 
Bowman  and  Seery  (1968) 


Reaction 

No. 

Mechanism 
1  or  2 

c2h2  +  c2h2 

t 

C4H2  +  h  +  H 

[4-1] 

1.2 

c2h2  +  0 

* 

ch2  +  CO 

[4-2] 

2 

C2h2  +  0 

t 

C2H  +  OH 

[4- 2b] 

1 

C2H2  +  OH 

t 

C2H  +  H20 

[4-3] 

1.2 

C2H2  +  H 

4- 

c2h  +  H2 

[4-4] 

1,2 

c2h  ♦  o2 

-* 

4- 

CO  +  CO  +  H 

[4-5] 

1,2 

C2H  +  0 

-¥ 

4- 

CO  +  CH 

[4-6] 

1 

CH2  +  02 

•¥ 

4- 

HCO  +  OH 

[4-7] 

2 

HCO  +  OH 

t 

CO  +  h2o 

[4-8] 

2 

CO  +  OH 

4- 

C02  +  H 

[4-9] 

1,2 

H  +  02 

* 

OH  +  0 

[4-10] 

1,2 

o  +  h2 

4- 

OH  +  H 

[4-11] 

1.2 

o  +  h2o 

4- 

OH  +  OH 

[4-12] 

1.2 

H  +  H20 

•+■ 

H2  +  OH 

[4-13] 

1,2 

H  +  OH  +  M 

-* 

4- 

H20  +  M 

[4-14] 

1,2 

H  +  H  +  M 

■+■ 

4- 

h2  +  h 

[4-15] 

1,2 
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0+H2  J  OH+H  [4] 

were  assumed  to  attain  equilibrium  (after  Schneider,  I960;  Hottel  et  al.^ 
1965).  Also  the  three-body  radical  recomb i nation  reactions: 


H  +  H  +  M  t  H2  +  M  [6] 

0  +  0  +  M  t  Og  +  M  [7] 

0  >  H  +  M  t  OH  +  M  [8] 

OH  +  H  +  M  *  H20  +  M  [9] 


were  allotted  to  proceed  In  the  forward  direction  only.  As  noted  pre¬ 
viously,  the  results  for  C^Hg  obtained  using  this  mechanism  were  not 
good. 

Chinltz  and  Baurer  (1966)  investigated  the  combustion  of  propane/ 
air  and  ethane/air.  They  postulated  a  combined  mechanism  containing 
69  basic  reaction  pairs  and  31  chemical  species;  the  ethane/air  mechanism 
appeared  as  a  subset  of  the  propane/air  mechanism.  The  complete  mech¬ 
anism  is  given  in  Table  5.  In  the  table  Reactions  [5-l]-[5-12]  pertain 
only  to  the  propane/air  combustion.  Using  an  analytical  definition  for 
ignition  delay  time  as  5%  of  the  final  tenperature  rise,  the  computed 
ethane  ignition  delay  times  were  fair  agreement  with  data  determined  in 
subsonic  diffusion  flames  using  vitiated  air.  However,  the  propane 
ignition  delay  times  compared  somewhat  less  well  with  some  data  obtained 
in  a  shock  tube.  The  experimental  ignition  delay  times  were  defined 
as  the  time  to  onset  of  luminosity,  and  thus  the  difference  in  definitions 
could  account  for  some  of  the  disagreement. 

For  this  combined  mechanism  our  division  of  the  combustion  into 
three  processes  fails.  In  particular.  Reactions  [5-18],  [5-19],  [5-21], 
[5-34]_[5_36],  and  [5-49]  all  produce  carbon  dioxide  rather  than  carbon 
monoxide  thus  invalidating  our  separation  of  partial  oxidation  from 
carbon  monoxide  destruction.  However,  Reaction  [5-61]  is  cited  by  Chinitz 
and  Baurer  as  the  predominant  carbon  monoxide  destruction  reaction, 
and  the  free  radical  reactions  are  essentially  the  same  as  in  the  pre¬ 
vious  mechanisms.  To  elucidate  the  coupling  between  the  partial  oxi¬ 
dation  and  carbon  dioxide  production  processes,  let  us  examine  the 
temperature-time  and  carbon  dioxide-time  histories  (cf.  Figs.  1  and  2) 
from  the  work  of  Chinitz  and  Baurer.  These  histories  show  a  brief 
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Table  5. 

The  combined  propane/air;  ethane/air 
combustion  mechanism  of  Chinitz  and  Baurer  (1966) 


Reaction 


C3H8  ♦  OH 

t 

^3^7  4  ^0 

[5-1] 

C3H8  +  0 

t 

C2Hg  +  HCHO 

[5-2] 

C3H8  ♦  H 

t 

C3H7  4  H2 

[5-3] 

C3H7  4  0 

t 

c2h5c::o  +  h 

[5-4] 

c3H7  +  0 

t 

C3H6  +  OH 

[5-5] 

C3H7  +  OH 

t 

C3H6  +  H20 

[5-6] 

C3H7  +  °2 

t 

C2H5CH0  +  OH 

[5-7] 

c3h6  +  o2 

-*■ 

4- 

CH3CH0  +  HCHO 

[5-8] 

C3H6  +  OH 

4- 

CH3CHO  +  CH3 

[5-9] 

c3h6  +  0 

% 

C3H4  ♦  h20 

[5-10] 

c3h6  +  0 

-*■ 

4- 

HCHO  +  C,H4 

[5-11] 

c3h4  ♦  o2 

4- 

CH3C0  CHO 

[5-12] 

c2h6  +  o 

4- 

C2H5  +  OH 

[5-13] 

C2H6  +  OH 

-*► 

4- 

c2h5  +  h2o 

[5-14] 

C2H6  -  H 

-*• 

4- 

c2h5  +  H, 

[5-15] 

C2H5CH0  +  OH 

4- 

c2h5co  +  H20 

[5-16] 

C2H5CH0  +  H 

t 

C2H5CO  +  H2 

[5-17] 

C2HrCH0  +  02 

4 

4- 

C2M5  4  C02  +  0H 

[5-18] 

c2h5co  +  0 

-4 

4- 

C2H5  +  C02 

[5-19] 
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Table  5.  Cont. 

Reaction 

No. 

C2H5C0  +  OH 

-4- 

C2H50H  +  CO 

[5-20] 

c2h5co  +  o2 

t 

C2H50  +  C02 

[5-21] 

C2H5C0  +  H 

t 

C2H5  +  CHO 

[5-22] 

c2h5oh  +  c2 

-*• 

CH30H  +  CO  +  H20 

[5-23] 

C2H5OH  +  H 

-4- 

C2H5  +  »2° 

[5-24] 

C2H5°  +  °2 

t 

CH30  +  CO  +  H20 

[5-25] 

C2H5°  +  °2 

-*• 

-4- 

CH30H  +  CO  +  OH 

[5-26] 

C2H50  +  H 

•4- 

C2H5  +  OH 

[5-27] 

^2H5  +  ®2 

*♦ 

-4- 

CH3CH0  +  OH 

[5-28] 

c2H5  +  0 

* 

C2H4  +  OH 

[5-29] 

C2H5  +  0 

■  ■> 

•4- 

CH3CH0  +  H 

[5-30] 

C2H5  +  OH 

-► 

■4- 

c2h4  +  h2o 

[5-31] 

CH^CHO  +  OH 

•4- 

ch3co  +  H20 

[5-32] 

CH3CH0  +  H 

HK 

4- 

ch3co  +  H2 

[5-33] 

ch3cho  +  o2 

CH3  +  C02  +  OH 

[5-34] 

ch3co  +  0 

t 

CH3  +  C02 

[5-35] 

ch3co  +  o2 

t 

CH30  +  C02 

[5-36] 

CH3CO  +  OH 

t 

CH30H  +  CO 

[5-37] 

CH3CO  +  H 

•* 

■4- 

CH3  +  CHO 

[5-33] 

CH30H  +  02 

-► 

■4- 

CO  +  H20  +■  H20 

[5-39] 

CH30H  +  H 

* 

ch3  +  H20 

[5-40] 

ch3o  +  02 

■4- 

CO  +  OH  +  H20 

[5-41] 

CH30  +  H 

-*• 

4- 

CH3  +  OH 

[5-42] 

C2H4  +  0 

-► 

4- 

CH3  +  CHO 

[5-43] 
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Table  5.  Cont. 


Reaction 


No. 


C2H4  +  0  t  C2H2  +  H20 

[5-44] 

C2H4  +  0H  t  CH3  +  HCHO 

[5-45] 

C2H2  +0  *  CgH  +  OH 

[5-46] 

C2H2  +  OH  t  C2H  +  H20 

[5-47] 

C2H  +  02  t  CO  +  CO  +  H 

[5-48] 

c2h  +  o2  ;  ch  +  co2 

[5-49] 

CH  +  02  t  CO  +  OH 

[5-50] 

CH4  +  OH  *  CH3  +  H20 

[5-51] 

CH4  +  0  t  CH3  +  OH 

[5-52] 

CH4  +  H  t  CH3  +  H2 

[5-53] 

CH3  +  0  t  HCHO  +  H 

[5-54] 

CH3  +  02  *  HCHO  +  OH 

[5-55] 

HCHO  +0H  J  CHO  +  H20 

[5-56] 

HCHO  +  H  +  CHO  +  H2 

[5-57] 

HCHO  +0  2  CHO  +  OH 

[5-58] 

CHO  +  0  t  CO  +  OH 

[5-59] 

CHO  +  OH  t  CO  +  H20 

[5-60] 

CO  +  OH  t  C02  +  H 

[5-61] 

OH  +  H2  t  H  +  H20 

[5-62] 

OH  +  OH  Z  H20  +  0 

[5-63] 

0  +  H2  t  OH  +  H 

[5-64] 

H  +  02  Z  OH  +  0 

[5-65] 

0  +  H  +  M  Z  OH  +  M 

[5-66] 

0+0+M  Z  02+H 

[5-67] 

H  +  H  +  M  Z  H2  +  H 

[5-68] 

H  +  OH  +  M  *  H20  +  H 


[5-69] 


4000 


The  effect  of  1 
(Chinitz  and  Ba 


16 


carbon  dioxide  and  temperature  overshoot  with  a  duration  of  approxi¬ 
mately  l(f4  seconds.  Chinitz  and  Baurer  attribute  these  overshoots 
to  the  Inability  of  the  CO-COg  equilibrium  (Reaction  [5-61])  to  match 
the  Initial  rapid  production  of  carbon  dioxide  by  the  reactions  men¬ 
tioned  above.  They  note  that  the  magnitude  and  existence  of  the  over¬ 
shoots  depend  on  the  values  of  some  of  the  less  accurately  known  rate 
constants.  If  indeed  the  overshoot  is  a  real  effect  (experimental 
verification  is  lacking)  its  duration  Is  short  and  its  effect  is  quickly 
damped.  In  fact,  the  overshoot  can  be  thought  of  as  indicating  the 
high  rate  of  the  partial  oxidation. 

The  common  feature  of  all  the  hydrocarbon  combustion  mechanisms 
presented  is  that  (with  the  exceptions  noted)  they  can  be  divided  Into 
three  basic  processes: 

1.  rapid  partial  oxidation  of  the  parent  hydrocarbon  to  carbon 
monoxide  and  water, 

2.  A  rate  limiting  step  of  carbon  monoxide  destruction  through 
Reaction  [24], 

3.  a  sequence  of  free  radical  production  reactions  not  Involving 
carbon. 

A  similar  division  was  made  by  Hammond  and  Hellor  (1970a, b)  and  by 
Edelman  (1970). 

In  the  present  work,  in  order  to  make  the  mechanism  applicable 
to  the  combustion  of  any  hydrocarbon  the  partial  oxidation  is  modelled 
using  tlte  following  general  reaction: 

CaHb  +  °2  -  aC0  +  |  H2°  C37] 

which  Is  assumed  to  have  an  infinitely  fast  forward  rate  compared  to 
the  rate  of  Reaction  [24],  The  necessity  of  making  such  a  gross 
characterization  of  the  process  becomes  apparent  when  one  considers 
the  extremely  complex  mechanism  which  was  developed  for  as  relatively 
simple  a  hydrocarbon  as  propane  (Chinitz  and  Baurer,  1966)  and  then 
envisions  the  complexity  of  the  mechanism  which  will  *•  'rely  be  required 
to  model  the  combustion  of  a  more  realistic  hydrocarbon  fuel  such  as 
kerosene.  Clearly,  our  present  knowledge  of  complex  hydrocarbon  com¬ 
bustion  is  insufficient  to  allow  development  of  accurate  mechanisms. 
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Edelman  (1970)  has  modelled  the  partial  oxidation  process  by  a  re¬ 
action  similar  to  Reaction  [37]  but  producing  carbon  monoxide  and 
diatomic  hydrogen  and  having  a  finite  but  empirically  determined 
rate.  A  similar  division  was  also  made  by  Marteney  (1970)  for  methane 
combustion.  In  agreement  with  our  postulate,  he  found  that  the  ''complete'' 
mechanism  results  agreed  with  those  using  an  infinitely  fast  partial 
oxidation. 

The  remainder  of  the  general  hydrocarbon  combustion  mechanism 
consists  of  the  single  carbon  monoxide  destruction  reaction  ([24])  and 
various  generally  accepted  free  radical  production  reactions  not  in¬ 
volving  carbon.  This  portion  of  the  mechanism  is  common  to  all  of  the 
cited  mechanisms  with  minor  deviations  in  the  free  radical  production 
reactions.  We  have  chosen  to  include  those  free  radical  reactions 
which  are  most  well  understood  and  have  accurate  accepted  rate  constants. 
Thus,  the  final  general  hydrocarbon  combustion  mechanism  is  shown  In 
Table  6. 

B.  The  nitric  oxide  formation  mechanism 

Since  the  nitric  oxide  emissions  from  gas  turbine  engines  ere  of 
current  Interest  from  a  pollution  viewpoint,  it  would  be  advantageous 
to  incorporate  a  nitrogen  oxidation  mechanism  into  the  hydrocarbon  com¬ 
bustion  mechanism.  Although  the  nitrogen  oxidation  chemistry  is  con¬ 
siderably  simpler  than  the  hydrocarbon  oxidation  chemistry,  only  recently 
has  a  well  defined  mechanism  been  developed. 

The  Zeldovich  mechanism  (Zeldovich  et  a!.,  1947)  consisting  of  the 
reactions: 


t 

NO  +  N 

[20] 

-4* 

4- 

NO  +  0 

[21] 

Is  generally  accepted  as  the  major  source  of  nitric  oxide  (Bowman,  1970; 
Newhall,  1969;  Camac  and  Feinberg,  1967;  Marteney,  1970;  Lavoie  et  al., 
1970;  Fletcher  and  Heywood,  1971).  The  only  point  of  disagreement  among 
these  various  authors  is  what  other  reactions  also  contribute  significantly 
to  the  formation  of  nitric  oxide.  This  point  must  be  examined  somewhat 
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more  closely. 

Bowman  (1970)  feels  that  only  Reactions  [20]  and  [21]  are  signifi 
cant.  Camac  and  Feinberg  (1967)  favor  including  the  reaction: 

N  +  0  +  K  t  NO  +  H  [22] 

for  high  temperatures  in  air  and  the  reaction: 

N2  +  02  t  NO  +  NO  [23] 

at  low  temperatures.  Marteney  (1970)  includes  only  Reactions  [20], 
[21],  and  [22]  in  his  analytical  work.  Newha.il  (1969)  initially  in¬ 
cludes  all  of  the  Reactions  [20]-[23]  and  two  additional  formation 
paths  involving  the  species  N02  and  N20  but  he  then  concludes  from 
his  result  that  the  Zeldovlch  mechanism  dominates.  Lavoie  et  al. 

(1970)  add  the  reaction: 

N  +  OH  t  NO  +  H  [38] 

giving  the  mechanism  shown  in  Table  7.  Fletcher  and  Heywood  (1971) 
also  use  the  mechanism  of  Table  7;  however,  both  investigators 
found  that  only  the  first  three  reactions  are  significant. 

The  basic  Zeldovich  mechanism  trust  clearly  be  incorporated  into 
any  reasonable  mechanism  of  nitric  oxide  formation.  BowTOn  (1970) 
established  from  his  shock  tube  studies  that  Reaction  [38]  could  be¬ 
come  significant  at  temperatures  below  2000°K,  but  we  have  not  included 
it  at  present.  In  deference  to  the  results  of  Camac  and  Feinberg 
(1967)  Reaction  [22]  is  included,  but  Reaction  [23]  is  excluded  since 
we  are  concerned  only  with  high  temperature  processes.  Our  final 
mechanism  thus  appears  in  Taole  8. 

Bowman's  (1970)  work  clearly  shews  an  advantage  which  our  total 
kinetic  model  gives  over  previous  models  which  assume  the  hydrocarbon 
chemistry  to  be  in  equilibrium  (Lavoie  et  al.,  1970;  Fletcher  and 
Heywood,  1971).  Note  the  discrepancy  resulting  from  the  equilibrium 
assumptions  (Fig.  3).  Bowman  accurately  attributes  this  deviation 
to  a  kinetic  oxygen  atom  overshoot  above  its  equilibrium  value. 
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Table  7 

Nitric  oxide  mechanism  of  Lavoie  et  al . 
(1970)  and  Fletcher  and  Heywood  (1971) 


Reaction 


N2  +  0 

t 

NO  +  N 

[20] 

N  +  02 

* 

NO  +  C 

[21] 

N  +  OH 

•4* 

4- 

NO  +  H 

[38] 

H  +  N20 

* 

N2  +  OH 

[39] 

o  +  n2o 

* 

n2  +  o2 

[40] 

o  +  n2o 

-4- 

NO  +  NO 

[41] 

Nitric  oxide  mechanism  (this  ref.) 


Reaction 


No. 


23 


C.  Hydrocarbon  combustion  rate  data 

The  determination  of  accurate  rate  constants  for  use  in  kinetic 
calculations  has  long  been  a  significant  problem.  For  the  reactions 
we  are  considering  the  comprehensive  surveys  of  Baulch  et  al.  (1968a, 
b,  1969a, b)  have  greatly  simplified  the  problem.  For  each  reaction 
pair  In  the  final  mechanism  we  have  taken  the  more  accurate  of  the 
forward  and  reverse  rate  constants  from  this  reference  where  possible 
and  expressed  it  In  the  modified  Arrhenius  form: 

kj  «  Aj  T6*  exp  (-Ej/RT)  (2-1) 

where  A^,  6^,  are  the  parameters  for  the  ith  elementary  reaction. 

The  less  accurate  of  the  rate  constant  pair  is  then  computed  from 
the  equilibrium  constant  (Stull,  1965)  and  the  other  rate  constant 
and  fitted  to  equation  (2-1).  The  values  of  the  forward  and  reverse 
rate  parameters  are  given  in  Table  9  for  all  reactions  In  the  combined 
hydrocarbon  combustion-nitrogen  oxidation  mechanism. 

D.  Formulation  of  the  governing  equations  for 
a  perfectly  stirred  reactor 

The  classical  concept  of  a  stirred  reactor  is  a  zone  of  uniform 
and  homogeneous  composition  and  temperature  through  which  a  steady 
flow  of  mass  passes  (Fig.  4).  Apy  mass  entering  this  zone  is  assumed 
to  be  instantaneously  mixed  throughout  the  volume  of  the  zone.  Thus 
inside  the  zone  mixing  Is  infinitely  fast,  and  kinetics  completely  con¬ 
trol  t:-'  reaction  rates.  Also,  since  the  zone  Is  perfectly  mixed 
the  composition  and  temperature  of  the  exit  flow  are  identical  to  the 
composition  and  temperature  Inside  the  zone.  With  this  basic  ground¬ 
work  we  can  proceed  to  formulate  the  governing  equations. 

Using  the  generalized  notation  of  Penner  (1957)  the  chosen  re¬ 
action  mechanism  may  be  represented  by: 

I  kfj  I 

*  xi  kt  Ev'iX  (j=l,2,...J) 

1*1  J  rj  1=1  1 


(2-2) 
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Table  9. 

Rate  parameters* 

Reaction  Forward  Reaction  Reverse  Reaction 


No. 

Afi 

6fi 

Efi 

Ari 

6ri 

Eri 

[1]* 

8.000  x  1014 

0. 

+45000 

1.071  x 

1012 

.39 

25751 

[2]+ 

2.190  x  1013 

0. 

+5150 

4.858  x 

1014 

-.20 

20727 

[3]+ 

2.240  x  1014 

0. 

+16800 

6.647  x 

1011 

.39 

-469 

[4]+ 

1.740  x  1013 

0. 

+19450 

7.737  x 

1012 

0. 

7467 

[5]+ 

5.750  x  1013 

0. 

18000 

1.169  x 

1012 

.20 

444 

r—i 

VO 

LJ 

5.006  X  1018 

-1.15 

0 

1.662  x 

1019 

-1.14 

103988 

[7]  + 

4.700  x  1015 

-.28 

0 

2.375  x 

1018 

.66 

119279 

[834 

5.300  x  1015 

0. 

-2780 

8.477  x 

1015 

0. 

99246 

C9]+ 

1.170  x  1017 

0. 

0 

8.629  x 

1018 

-.19 

119565 

[20]+ 

6.192  x  1013 

.10 

75241 

3.100  x 

1013 

0. 

334 

[21  f 

6.430  x  109 

1.00 

6250 

3.661  x 

108 

1.16 

37847 

[22  f 

14** 

6.450  x  10'^ 

-.50 

0 

1.856  x 

1016 

-.72 

150876 

[24]+ 

5.600  x  1011 

0. 

1080 

1.455  x 

1018 

-1.19 

27C47 

*  uni ts , 

g-mole,  cc,  °K, 

,  cal 

f  forward  parameters  from  Jenkins  et  al,  (1967) 

+  forward  parameters  from  Baulch  et  al .  (1968a,b,  1969a, b) 

**  Due  to  a  coding  error  this  valve  is  lower  by  a  factor  of  100  than 
that  of  Baulch,  et  al.;  however,  calculations  reveal  that  this 
error  has  a  negligible  effect  ori  the  predictions  to  be  presented 
in  Chapter  V. 
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where  denotes  the  i**1  species,  v!j  its  stoichiometric  coefficient 

as  a  reactant,  and  vjj  its  stoichiometric  coefficient  as  a  product, 

where  there  are  I  species  and  J  reactions  (here  I  =  12,  J  =  14).  The 

symbols  k~.  and  k  .  denote  the  forward  and  reverse  rate  constants  for 
fj  rj 

the  j*-n  reaction,  respectively. 

The  species  conservation  equations  are  given  in  general  by: 

mtl  -  mi2  *  -  mj"  V  (i=l,2,...I)  (2-3) 

where  m.  is  the  total  mass  flow  of  the  1^  species,  irJ"  is  the  volum- 

I  1.L  ' 

etric  mass  production  rate  of  the  iin  species  inside  the  reactor  due 
to  chemical  reaction,  and  V  is  the  reactor  volume.  The  subscripts  1 
and  2  denote  inlet  and  exit  (also  internal)  conditions,  respectively. 
After  Jones  and  Prothero  (1958)  equation  (2-3)  can  be  recast  as: 

m1  *  * 

y  (cr^i  -  =  “  yj —  (1=1,2, ...I)  (2-4) 

X  L 

where  a,  is  the  concentration  of  the  itnspecies  (mole  i/gm),  m  is  the 
total  mass  flow  rate,  and  is  the  molecular  weight  of  the  1  species. 

Assuming  local  thermal  equilibrium,  a  mixture  of  thermally  perfect 
gases,  and  constant  reactor  pressure,  the  source  term  in  equation  (2-4) 
can  be  written  explicitly  giving: 


V  [ail  "  °i2^ 


jf1  (vij  "  vij)  tkfj  i=1  ai2)Vij 

1  X," 

-  k  .  tt  (p  oi?)vij]  (i»l,2,...I)  (2-5) 

rj  i=1  12 


and  the  gas  law: 


P 


P 

- j - 

RT0  E  0., 
1  1=1  U 


(2-6) 


Assuming  negligible  changes  in  kinetic  and  potential  energy  and 
no  shaft  work,  the  energy  equation  can  be  written  as: 
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I 

i^1<CT12h12  "  °i!hi^  “  q  (2‘7) 


where  h^  Is  the  total  (chemical  +  sensible)  enthalpy  of  the  1th  species 
(cal/mole  1)  and  q  Is  the  heat  added  (cal/g  of  flow). 

The  set  of  equations  (2-5),  (2-6),  and  (2-7)  number  I  +  2,  and 
there  are  1  +  2  dependent  variables  (a.j2,  ^2’  therefore  the  problem 
is  well  posed.  A  computer  program  was  written  to  solve  the  set  of 
equations  simultaneously  (cf.  Appendix  A).  We  will  now  present  a  dis¬ 
cussion  of  the  predictions  obtained  from  the  program. 

E.  Comparison  of  the  predicted  results  with 
experimental  data 

The  true  measure  of  the  usefulness  of  this  perfectly  stirred  re¬ 
actor  analysis  in  building  an  overall  gas  turbine  combustor  model  will 
be  Its  ability  to  predict  the  combustion  performance  of  physical  systems. 
Any  error  Incurred  in  predicting  the  performance  of  a  single  perfectly 
stirred  reactor  will  certainly  manifest  itself  in  an  error  In  the  over¬ 
all  predictions.  For  this  reason  we  have  checked  the  accuracy  of  our 
analysis  In  two  ways:  the  first  is  an  asymtotic  comparison  with  equil¬ 
ibrium  results,  and  the  second  is  a  comparison  with  experimental  data 
for  a  single-stage  reactor. 

It  Is  expected  that  the  kinetic  predictions  approach  equilibrium 
values  as  the  mass  flow  rate  Is  decreased  or  the  inlet  tenperature  in¬ 
creased,  To  check  this  correspondence  Fig.  5  was  prepared.  In  this 
and  subsequent  figures  completeness  of  combustion  e  is  defined  by: 


c 


XC0?,2 

xC02,2  +  XC0,2 


(2-8) 


where  x^2  is  the  reactor  and  exit  mole  fraction  of  the  indicated  species. 
The  curve  denoted  as  equilibrium  Is  the  result  of  an  independent  com¬ 
putation  for  the  same  conditions  based  on  accepted  values  of  the  equil¬ 
ibrium  constants  (Stull,  1965).  The  equilibrium  computer  program  is 
described  In  Appendix  B.  In  the  figure  the  kinetic  results  are  seen 
to  tend  towards  equilibrium  as  the  flow  rate  per  unit  volume  is  decreased 
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(ror  a  constant  pressure  of  one  atm).  The  results  also  exhibit  the 
correct  inlet  temperature  trend,  i.e.  for  constant  flow  rate  per  unit 
volume  and  constant  pressure  the  results  tend  toward  equilibrium  values 
with  increasing  inlet  temperature.  Finally,  at  a  given  loading  the 
completeness  of  combustion  e  (a  measure  of  combustion  efficiency)  de¬ 
creases  with  Increasing  inlet  temperature.  This  trend,  observed  in 
our  earlier  results  (Hammond  and  Mellor,  1970a),  is  consistent  with 
the  explanation  of  increased  carbon  dioxide  dissociation  at  higher  re¬ 
actor  temperatures. 

In  order  to  further  define  the  accuracy  of  our  single  stage  stirred 
reactor  program  for  various  fuels,  comparisons  were  made  with  experi¬ 
mental  data  for  the  combustion  of  carbon  monoxide/air  In  a  spherical 
combustor  (Nerheim,  1960),  methane/air  in  conical  and  cylindrical  com¬ 
bustors  (Morgan,  1967),  and  propane/air  In  a  spherical  combustor 
(Schneider,  1960).  These  coiparisons  are  shown  in  Figs.  6-11. 

Considering  first  the  carbon  monoxide/air  data  (Figs.  6-9),  one 
would  expect  the  theory  to  be  most  accurate  for  this  fuel  since  our 
hydrocarbon  combustion  scheme  is  based  on  a  rate  limiting  carbon  mon¬ 
oxide  combustion  mechanism.  This  comparison  should  thus  provide  a  test 
of  the  validity  of  this  generally  accepted  carbon  monoxide  combustion 
mechanism.  Fig.  6  reveals  a  rather  large  discrepancy  between  the  pre¬ 
dicted  and  experimental  values,  i.e.  theory  predicts  blow-out  while 
the  experimental  values  clearly  indicate  the  presence  of  combustion. 

The  discrepancy  In  Fig.  6  results  from  our  use  of  Reaction  [24]  as  the 
only  carbon  monoxide  destruction  reaction.  In  the  presence  of  high 
water  concentrations  (and  thus  an  ample  supply  of  hydroxyl  radicals) 
Reaction  [24]  is  truly  dominant,  but  as  the  amount  of  water  is  reduced 
other  reactions  have  been  shown  to  become  significant  (Newhall,  1969). 
Therefore,  for  "dry"  carbon  monoxide/air  combustion  the  following  re¬ 
actions  should  be  considered: 

co  +  o2  t  co2  +  0 

CO  +  0  +  M  t  C02  +  M 

Since  these  reactions  are  not  included  in  our  mechanism  the  rate  of 
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FIGURE  7  THEORETICAL  STIRRED  REACTOR  PERFORMANCE 
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FIGURE  8  THEORETICRL  STIRRED  REACTOR  PERFORMANCE 
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carbon  monoxide  destruction  will  be  underestimated  as  the  hydroxyl 
radical  concentration  decreases.  Fortunately,  for  the  hydrocarbon/ 
air  systems  we  seek  to  model  sufficient  hydroxyl  radical  concentrations 
are  insured  by  the  hydrogen  content  of  the  fuel.  Thus,  the  discrepancy 
evidenced  is  considered  to  result  primarily  from  a  decrease  in  the  rate 
due  to  neglecting  Reactions  [25]  and  [26]  and  a  reduction  in  the  rate 
of  Reaction  [24]  due  to  a  lack  of  hydroxyl  radicals  and  will  therefore 
not  be  important  for  hydrocarbons. 

Mixing  in  the  physical  system  would  play  some  role  in  the  perfor¬ 
mance  degradation  of  the  experimental  data  of  Figs.  8-9,  but  the  large 
scatter  in  the  experimental  data  makes  definite  conclusions  difficult. 

The  methane/air  combustion  data  are  presented  in  Figs.  10-11. 

Although  the  discrepancy  between  theory  and  experiment  still  persists, 
its  magnitude  is  somewhat  smaller  than  for  the  previous  data.  In  both 
cases  the  trends  predicted  are  completely  correct.  The  predicted  per¬ 
formance  for  CH4  exceeds  the  experimental  performance  as  a  result  of 
assuming  an  infinitely  fast  partial  oxidation  step.  It  has  been  reason¬ 
ably  well  established  that  Reaction  [2-1]  constitutes  a  second  rate 
limiting  step  in  the  methane  combustion  process  (Morgan,  1367).  Therefore, 
the  partial  oxidation  process  cannot  have  an  infinite  rate,  and  as  a 
result  the  predicted  values  are  high.  Mixing  limitations  as  always 
could  further  reduce  the  experimental  performance. 

In  general  Reaction  [2-1]  has  not  been  found  to  be  rate-limiting 
in  the  combustion  of  hydrocarbons  higher  than  methane.  The  assumption 
of  an  infinitely  fast  partial  oxidation  should  thus  improve  for  higher 
hydrocarbons.  The  data  for  propane/air  combustion  shown  in  Fig.  11 
appear  to  substantiate  this  hypothesis.  The  agreement  is  extremely 
good  over  the  range  of  the  data,  except  for  a  slight  deviation  in  the 
region  of  low  loadings.  This  deviation  is  most  likely  a  result  of  in¬ 
creased  segregation  in  the  physical  system  as  the  mixing  rate  decreases 
with  decreasing  flow  rate. 

To  establish  the  ability  of  the  kinetic  analysis  for  a  single 
stage  perfectly  stirred  reactor  to  realistically  predict  nitric  oxide 
fornatlon,  its  predictions  were  compared  with  the  aforementioned  equil¬ 
ibrium  calculations  (Fig.  12).  In  all  cases  the  predicted  nitric  oxide 


concentration  increases  with  increasing  inlet  temperature,  and  as  the 
flow  rate  is  decreased  the  results  tend  toward  equilibrium.  At  very 
low  flows  the  equilibrium  and  kinetic  curves  essentially  coincide*,  in 
fact,  the  discrepancy  between  them  is  a  measure  of  the  numerical  accuracy, 
and  is  less  than  12. 

In  view  of  the  good  correspondence  between  theory  and  experiment 
for  propane  and  the  fair  correspondence  for  methane,  the  analysis  will 
be  used  "as  1s“  for  hydrocarbons  above  and  Including  propane.  Investi¬ 
gation  of  a  more  appropriate  model  for  partial  oxidation  is  continuing; 
however,  a  vast  number  of  reactions  which  would  seriously  increase  com¬ 
putational  cost  Is  clearly  undersirable. 

In  all  of  the  comparisons  above  mixing  has  been  blamed  for  at  least 
a  portion  of  the  discrepancy  observed  between  theory  and  experiment. 

To  further  define  the  exact  effects  of  mixing  we  now  will  investigate 
It  in  more  detai  1 . 


CHAPTER  III 


THE  INFLUENCE  OF  MIXING  ON  THE  COMBUSTION  PROCESS 

The  performance  of  any  physical  combustion  system  will  be  deter¬ 
mined  by  the  rates  of  two  processes:  chemical  reaction  and  mixing. 

In  the  previous  chapter  systems  having  infinite  mixing  rates  were  con¬ 
sidered.  In  this  chapter  we  turn  to  an  investigation  of  the  effects 
of  a  finite  mixing  rate  on  such  systems  both  through  Its  effect  on  the 
residence  time  distribution  and  through  segregation.  Considering  a 
zone  of  given  volune,  the  residence  time  distribution  gives  the  time 
which  each  fluid  element  (which  may  have  finite  or  molecular  size)  has 
been  In  the  system  and  the  degree  of  segregation  Indicates  the  freedom 
with  which  these  fluid  elements  may  intermix.  Each  of  these  concepts 
will  be  considered  in  turn,  but  first  It  Is  advantageous  to  consider 
the  turbulent  mixing  process  In  general. 

In  the  theory  of  micro-volume  combustion  reviewed  by  Hamnond  and 
Mellor  (1970a, b)  It  was  stated  that  the  turbulent  mixing  process  could 
be  characterized  by  the  formation  and  subsequent  decay  of  eddies  in 
the  flow,  and  this  concept  Is  useful  here.  These  eddies  are  visualized 
as  reasonably  discrete  aggregates  of  fluid  which  retain  their  Identity 
for  some  lifetime  determined  by  the  turbulent  parameters  of  the  flow. 

If  the  properties  Inside  each  Individual  eddy  are  assumed  to  the  un¬ 
iform  (in  the  limit  of  zero  eddy  volume  the  assumption  Is  exact),  then 
the  segregation  of  the  fluid  is  directly  related  to  the  scale  of  the 
eddies,  and  the  decay  of  segregation  is  directly  related  to  their 
lifetimes. 

If  In  a  certain  flow  the  eddies  posses  Infinite  lifetimes,  the  flow 
Is  said  to  be  completely  segregated,  and  the  fluid  is  said  to  behave  as  a 
macrofluid.  On  the  other  hand,  if  the  eddy  lifetime  is  zero,  the  flow  is 
said  to  have  no  segregation,  and  the  fluid  is  said  to  behave  as  a  micro¬ 
fluid.  The  terms  macromixing  and  micromixing  have  beer,  used  to  describe  the 
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accompanying  mixing  processes  (Levensplel,  1962).  Fig.  13  shows  a 
physical  Interpretation  of  a  microfluid  and  a  macrofluid.  A  quantita¬ 
tive  definition  of  segregation  was  given  by  Danckwerts  (1958b). 

Therefore,  several  items  must  be  known  to  completely  define  the 
mixing  performance:  the  residence  time  distribution  of  the  eddies,  the 
rate  of  eddy  intermixing,  and  the  initial  distribution  of  reactants 
within  the  various  eddies.  Firstly,  we  will  consider  the  determination 
of  the  residence  time  distribution,  secondly,  the  initial  distribution 
of  reactants,  and  finally,  the  rate  of  eddy  intermixing. 

A.  The  determination  of  residence  time  distributions 

The  residence  time  behavior  of  a  system  is  generally  described 
In  terms  of  four  distribution  functions.  These  descriptions  were  first 
presented  in  part  by  Gilliland  and  Mason  (1949,  1952)  and  expanded  by 
Danckwerts  (1953)  and  Levensplel  (1962).  Any  flow  reactor  is  character¬ 
ized  by  a  certain  mean  residence  time  defined  by: 

T  a  pV/m  (3-1) 

where  "p  is  the  mean  fluid  density  Inside  the  reactor,  V  the  reactor 
volume,  and  m  the  mass  flow  rate.  For  further  analysis  it  is  convenient 
to  nondime  ns  Iona  lize  time  on  the  mean  residence  time,  thus: 

0  =  t/x  (3-2) 

where  t  is  dimensional  time  and  0  Is  called  the  reduced  time.  Except 
in  certain  extreme  cases  t  is  not  known,  since  p  is  not  known  until  the 
Internal  gas  composition  is  known  from  the  solution;  therefore  a  distri¬ 
bution  of  residence  times  must  be  considered. 

Intuitively  one  would  not  expect  all  of  the  fluid  within  the  vessel 
to  possess  Identically  the  same  residence  time;  therefore  the  internal 
age  distribution  function  I  is  defined  such  that  IdO  is  the  fraction  of 
material  inside  the  reactor  wi th  residence  time  between  0  and  0  +  do. 
Since  each  fluid  element  must  possess  a  residence  time  between  zero  and 
Infinity,  Integral  of  I  over  all  time  must  be  unity. 

In  an  analogous  manner,  E,  the  exit  residence  time  distribution 


function.  Is  defined  such  that  Ed0  is  the  fraction  of  material  leaving 
the  reactor  with  residence  time  between  0  and  f>  +  do.  E  must  satisfy 
the  same  normalization  condition  as  I.  Note  that  in  the  case  of  com¬ 
plete  or  no  segregation  the  E  function  can  be  used  as  a  weighting  func¬ 
tion  to  integrate  the  set  of  equations  defining  the  chemical  reaction 
rate,  thereby  giving  the  overall  performance  (assuming  the  initial  re¬ 
actant  distribution  Is  known).  The  performance  of  partially  segregated 
systems  cannot  be  simllarlly  obtained,  since  the  reaction  rate  then 
also  depends  on  the  rate  of  eddy  intermixing. 

The  E  and  I  distribution  functions  are  related  to  two  other  func¬ 
tions  which  can  be  evaluated  experimentally.  Consider  now  some  tracer 
fluid  whose  presence  can  be  easily  detected.  The  F  function  is  de¬ 
fined  as  the  nondlmensional  tracer  concentration  response  to  a  unit 
step  function  tracer  input.  Similarily,  the  C  function  is  the  non- 
dimensional  tracer  concentration  response  to  a  unit  delta  function 
tracer  input.  The  derivation  of  the  following  relations  among  these 
functions  Is  straightforward  and  will  not  be  repeated  here  (for 
details  see  Levenspiel,  1962  and  Danckwerts,  1953); 


F  +  I 


(3-3) 


E  = 


_  dF 
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(3-4) 

(3-5) 


The  1,  F,  and  E  functions  for  the  two  idealized  cases  of  per¬ 
fectly  stirred  flow  and  plug  flow  are  shown  in  Figs.  14  and  15,  re¬ 
spectively,  both  for  a  unit  step  input  function.  As  pointed  out  by 
Levenspiel  (1962),  and  E  and  F  functions  are  not  sufficient  to  as¬ 
certain  the  segregation  of  the  system.  Only  in  the  extreme  of  com¬ 
plete  segregation  can  the  residence  time  distribution  derived  with  the 
aid  of  the  F  and  C  functions  be  interpreted  as  the  eddy  residence 
time  distribution,  because  only  in  this  case  does  an  eddy  retain  its 
identity  as  a  discrete  fluid  element. 

Two  basic  models  have  been  developed  to  treat  the  effects  of 


Fig.  14  I,  F,  and  E  functions  for  a  perfectly  stirred 
reactor  (after  tevenspie’,  196?) 
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Imperfect  mixing  on  the  residence  time  distribution  functions.  The 
first  is  called  the  dispersion  model  (Levenspiel,  1962)  and  treats 
mixing  as  a  deviation  from  ideal  plug  flow.  Considering  a  reactor 
of  characteristic  length  l  through  which  fluid  passes  with  a  bulk  vel- 
oc  ■  U,  the  action  of  turbulent  and  molecular  transport  wi 11  certainly 
act  to  produce  some  mixing  in  the  flow  direction  of  fluid  elements 
having  different  residence  times.  Oanckwerts  (1953)  derived  an  analytical 
form  for  the  F  function  from  which  the  E  function  can  be  obtained  as: 


E(0) 


1  +  0 


r  T 


(3-6) 


where  D  is  a  turbulent  mass  diffusivity.  Eqn.  3-6  applies  for  the  case- 
4D/LU  <  <  1.  The  form  of  the  E  function  Is  more  clearly  shown  in  Fig. 

16,  and  the  predicted  deviation  from  plug  flow  can  be  seen.  The  figure 
shows  that  for  D/LU  =  0  the  E  function  is  that  of  plug  flow,  and 
for  Increasing  D/LU  the  function  deviates  more  strongly  from  that  for 
plug  flow.  The  complete  range  of  mixing  from  plug  to  perfectly  stirred 
flow  cannot  be  predicted  because  the  solution  is  restricted  to  4D/LU  <  <  1 
as  noted. 

The  second  model,  called  the  tanks-in-series  model,  approximates 
mixing  as  a  deviation  from  perfectly  mixed  flow.  After  Levenspiel  (1962), 
the  reactor  to  be  modelled  is  divided  into  a  series  of  equal  volume 
zones  which  are  perfectly  stirred  reactors.  The  predicted  E  function  is: 

E(O)  =  fR~T5Texp  (3'7) 

where  N  is  the  number  of  zones.  The  E  function  for  this  model  is  shown 
in  Fig.  17.  Note  that  the  full  range  of  behavior  from  perfectly  stirred 
(N«l)  to  plug  flow  (U+“)  can  be  predicted.  Therefore,  this  model  is 
better  from- a  practical  viewpoint  than  the  dispersion  model. 

The  previous  two  models  are  used  in  practice  either  by  arbitrarily 
setting  the  value  of  D/LU  and  N  or  by  using  these  values  as  paraneters 
to  curve  fit  experimental  data.  In  general,  these  models  alone  are  not 


16  The  E  functions  predicted  by  the  dispersi 
(after  Levenspiel,  1962). 
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Fig.  17  Tha  E  functions  predicted  by  the  tanks-i n-serles 
mode  1 . 


sufficient  to  model  a  coirfeustor.  By  using  various  combinations  of  these 
models  in  conjunction  with  pure  plug  and  perfectly  stirred  flow,  bypass 
flow,  and  recirculation  flow  one  can  approximate  a  known  residence  time 
distribution.  Two  major  difficulties  remain,  however:  1)  these  re¬ 
sults  are  useful  only  if  complete  or  m  segregation  Is  assumed;  2)  un¬ 
less  the  volumes  and  flows  are  arbitrarily  assigned  a  model  cannot  be 
formulated  without  testing  a  physical  system. 

As  a  direct  result  of  the  two  difficulties  mentioned  above,  all 
current  models  attempting  to  predict  the  various  residence  time  dis¬ 
tribution  functions  for  practical  systems  are  either  purely  postulative 
or  of  a  semi -empirical  nature.  We  propose  to  discuss  only  the  most 
promising  models  here. 

Beer  and  Lee  (1965)  derived  an  empirical  residence  time  distribu¬ 
tion  for  the  IJmuiden  tunnel  furnace  using  a  1/10  scale  waxer  model 
.iith  salt  tracer.  The  system  was  essentially  a  confined,  cylindrical, 
swirling  jet.  They  were  able  to  approximate  their  experimental  F 
functions  by  a  perfectly  stirred  and  a  plug  flow  reactor  in  series. 

The  fraction  of  total  residence  time  (and  hence  total  volume)  devoted 
to  perfectly  stirred  flow  was  determined  as  a  function  of  the  swirl 
number: 


where  Gx  is  the  total  swirl  angular  momentum,  Gv  is  the  total  axial 

<J)  A 

linear  momentum,  and  rQ  is  the  inlet  jet  radius.  This  function  is 
shown  in  Fig.  18.  With  no  sv/irl  the  ratio  has  a  value  of  0.5  end  shews 
a  rapid  decrease  with  increasing  S  up  to  a  value  of  S  =  1.0.  For 
S  >  1.0  the  fraction  devoted  to  perfectly  stirred  flow  increases 
steadily.  Beer  and  Lee  explain  their  curve  from  the  following  sequence 
of  physical  events  involving  a  wall  and  a  central  core  recirculation 
zone.  Initially  a  recirculation  zone  is  formed  along  the  chamber  walls 
due  to  the  entrainment  effect  of  the  jet.  Viewing  the  jet  as  devoted 
primarily  to  plug  flow  and  the  recirculation  regions  to  perfectly 
stirred  flow,  the  initial  increase  in  swirl  serves  only  to  increase 
the  jet  spread  angle  reducing  the  wall  recirculation  zone  volume  with- 


s 

Fig.  18  Residence  time  fraction  devoted  to  perfectly  stirred 
flow  (Beer  and  Lee,  1965). 
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out  forming  the  central  recirculation  zone.  Hence,  a  drop  In  per¬ 
fectly  stirred  volume  occurs.  At  S  :  1  the  central  recirculation 
zone  starts  to  form  (Chlgler  and  Chervlnsky,  1967)  due  to  the  radial 
pressure  gradient  and  Increases  In  size  with  Increasing  $;  therefore 
the  net  perfectly  stirred  volume  Increases.  Using  this  model  Beer 
and  Lee  obtained  reasonable  agreement  with  experiment. 

Pratt  (1970)  also  considered  a  confined,  swirling  jet.  He  used 
the  empirical  recirculation  model  of  Thrlng  and  Newby  (1953)  as  mod¬ 
ified  by  Field  et  al .  (1965)  and  the  mass  entrainment  function  of 
Kerr  and  Frazer  (1965).  By  considering  differential  equations  for 
the  conservation  of  a  tracer  substance  he  Indicated  how  to  evaluate 
the  F  function  for  the  jet.  He  was  able  to  predict  a  characteristic 
delay  time  8d  which  is  the  time  of  first  change  in  the  F  function  to 
some  input  and  noted  that  It  should  correspond  to  the  volume  fraction 
devoted  to  plug  flow  as  defined  by  Beer  and  Lee  (1965).  His  work  was 
restricted  to  low  swirl  (S  <  1),  but  fair  agreement  with  the  data  of 
Beer  and  Lee  was  obtained  In  that  range. 

Crowe  and  Pratt  (1970)  studied  the  conical  reactor  of  Jain  and 
Spalding  shown  In  Fig.  19.  Following  Jain  and  Spalding,  they  defined 
a  certain  fraction  R  of  the  Inlet  flow  which  recirculates  Inside  the 
reactor.  Using  the  analysis  of  Pratt  (1970)  they  predicted  the  F 
functions  of  the  system  (Fig.  20).  The  net  result  Is  a  one  parameter 
model  for  the  residence  time  distribution  since  R  cannot  be  predicted. 
The  advantage  of  this  model  is  that  it  predicts  a  characteristic 
delay  time . 

The  two  oarameter  recirculation  model  of  Hottel  et  al.  (1958) 
was  discussed  in  previous  work  (Hammond  and  Mellor,  1970a,  b).  It 
was  shown  to  predict  a  wide  variation  in  mean  residence  time  but  al¬ 
ways  predicts  a  zero  characteristic  delay  time. 

In  simmary,  we  see  that  the  residence  time  distribution  functions 
can  be  approximated  using  a  number  of  various  models  having  one  or 
more  parameters,  but  that  no  priori  analyses  exist.  We  now  consider 
determining  the  Initial  distribution  of  reactants. 
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B.  The  determination  cf  the  initial 
distribution  of  reactants 

Currently  only  two  attenpts  have  been  made  to  model  the  initial 
distribution  of  reactants.  The  first  is  a  three  parameter  permutation 
of  the  recirculation  model  of  Hottel  el  al.  (1958)  which  has  an  add¬ 
itional  fuel  injection  point  (Fig.  21).  Since  each  reactor  is 
assumed  to  be  micrcmixed,  in  essence  two  large  "eddies"  are  considered. 

Fletcher  and  Heywocd  (1971)  modelled  the  primary  zone  of  a  gas 
turbine  combustor  as  a  single,  completely  segregated  perfectly  stirred 
reactor.  At  the  reactor  inlet  a  certain  normal  distribution  of  equiv¬ 
alence  ratio  was  assumed,  and  both  the  mean  and  standard  deviation 
could  be  varied  arbitrarily,  A  normal  distribution  was  chosen, apparently 
from  consideration  of  fuel  spray  characteristics. 

C.  The  determination  of  the 
rate  of  eddy  intermixing 

The  rate  of  eddy  intermixing  is  important  in  evaluating  reactions 
in  partially  mixed  flows  since  it  is  the  single  mechanism  for  bringing 
reactants  together.  The  eddy  mixing  or  lifetime  problem  has  been 
discussed  by  Hammond  and  Mel! or  (1970a,b). 

Fletcher  and  Heywood  (1971)  attacked  the  problem  by  assuming  that 
the  eddy  decays  0''er  some  arbitrary  length,  and  at  the  end  of  that  length 
the  flow  is  micromixed.  This  was  accomplished  by  forcing  the  standard 
deviation  of  their  species  concentration  distributions  to  zero. 

Summarizing  the  net  contributions  of  all  these  mixing  models  we 
see  that: 

1)  the  exit  residence  time  distribution  function  can  be  determined 
working  from  semi -empirical  equations; 

2)  the  initial  distribution  of  reactants  is  usually  assumed  to  be 
uniform,  but  in  some  cases  an  educated  guess  at  a  distribution  can 
be  made; 

3)  some  sort  of  eddy  intermixing  function  must  be  assumed. 

He  must  now  consider  how  this  limited  knowledge  of  mixing  may  be 
applied  to  aid  performance  predictions  in  a  practical  combustion  system. 


CHAPTER  IV 

DEVELOPMENT  OF  AN  ANALYTICAL 
GAS  TURBINE  COMBUSTOR  MODEL 


The  previous  two  chapters  have  been  devoted  to  detailed  Investiga¬ 
tions  of  the  two  rate  processes  which  interact  to  define  the  overall 
behavior  of  any  physical  combustion  system.  It  now  remains  to  assemble 
the  results  of  these  two  studies  into  an  integrated  and  workable  form 
for  application  to  the  problem  of  analytically  modelling  gas  turbine 
combustors.  To  elucidate  the  most  realistic  merging  of  mixing  and 
kinetics,  the  physical  characteristics  of  gas  turbine  cont>ustors  will 
be  examined. 

A.  Pertinent  physical  and  flow  processes 
in  gas  turbine  combustors 

The  configuration  of  a  representative  gas  turbine  combustor  is 
shown  in  Fig.  22  and  the  flow  pattern  is  schematically  indicated.  For 
tubular  and  tuboannular  configurations  the  figure  shows  a  centerline 
cross-secticii  o.‘  the  axisymmetric  configuration;  for  an  annular  con¬ 
figuration  the  section  is  taken  on  a  plane  passing  through  the  engine 
centerline  and  the  combustor  is  considered  to  be  two-dimensional.  For 
either  case  the  essential  features  of  the  flow  pattern  and  the  discussion 
which  follows  are  identical. 

Consider  first  the  primary  zone  of  the  cortnistor  which  is  taken 
as  the  volume  contained  in  the  flame  tube  from  the  start  of  the  com¬ 
bustor  up  to  the  centerline  of  the  first  row  of  penetration  holes. 

There  may  or  may  not  be  some  addition  of  film  cooling  air  in  this  zone. 
Examining  the  schematic  flow  pattern  in  this  zone  and  the  experimental 
flow  pattern  of  Hiett  and  Powell  (1%2)  shown  in  Fig.  23,  it  is  apparent 
that  the  primary  zone  contains  a  complex  recirculating  flow  pattern. 

Also  all  of  the  fuel  is  injected  in  this  zone  and  must  mix  with  the 
available  air  and  hot  gases. 


primary  combustion  ,  secondary  combuation  .secondary  combustion  dilution  zone 
zone  I  zone  zone 
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r-  recirculation  zone 


Fig,  23  Experimental  primary  zone  flow  pattern 
(Hiett  and  Powell,  1962) 
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Heywood  et  al .  (1970)  view  the  fuel  evaporation  and  mixing  pro¬ 
cess  as  proceeding  with  negligible  combustion  until  the  fuel  is  in 
the  vapor  phase  entirely.  They  postulate  an  extremely  fuel  rich  re¬ 
gion  close  to  the  injector  where  the  atomized  fuel  spray  droplets 
evaporate.  No  combustion  is  allowed  in  this  region  due  to  a  lack  of 
air  as  well  as  the  postulate  that  the  fuel  drops  have  such  a  large 
velocity  relative  to  the  gas  phase  that  a  droplet  type  diffusion  flame 
cannot  be  sustained.  The  size  and  shape  of  the  vaporization  zone  will 
of  course  depend  on  the  injector  and  swirler  characteristics.  Both 
the  swirler  air  and  recirculating  hot  gases  entrain  vaporized  fuel  from 
this  region,  and  thus  combustion  takes  place  predominantly  in  the  hom¬ 
ogeneous  mode.  The  initial  distribution  of  vaporized  fuel  over  the 
combustor  cross  section  will  determine  the  initial  concentration  distri¬ 
bution  function.  Heywood  and  Fletcher  (1971)  assumed  this  to  be  a 
normal  distribution  as  mentioned  previously;  unfortunately  this  is  the 
only  available  solution  at  present.  In  consensus,  this  treatment  of 
vaporization  and  mixing  implies  that  the  actual  rate  of  vaporization 
is  not  important  and  that  the  only  significant  parameter  is  the  initial 
distribution.  The  effects  of  finite  evaporation  rates  have  been  studied 
qualitatively  by  Lefebvre  (1966,  1967). 

Although  the  actual  primary  zone  flow  pattern  is  heavily  dependent 
on  the  physical  configuration  of  the  flame  tube  in  the  area  (the 
swirler- injector  assembly  and  the  first  row  of  secondary  air  penetration 
holes),  three  distinct  recirculation  modes  can  usually  be  defined; 

1.  a  reverse  flow  zone  established  along  the  flame  tube  wall  as 
a  result  of  the  abrupt  increase  ir.  diameter  of  the  flame  tube 
over  that  of  the  swirier; 

2.  a  central  core  reverse  flow  zone  established  by  the  adverse 
centerline  pressure  gradient  set  up  by  swirl  (in  some  cases 
this  zone  becomes  annular  as  a  result  of  injector  interference); 

3.  a  centerline  reverse  flow  at  the  downstream  end  of  the  primary 
.  zone  resulting  from  the  central  impingement  of  the  first  row 

of  secondary  air  penetration  jets  (in  some  cases  this  merges 
with  the  reverse  flow  central  core , . 

A  similar  division  is  made  by  Way  (1955).  To  reasonably  model  the 
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primary  zone  each  of  these  flow  reversals  must  be  investigated.  \ 

The  first  type  of  recirculation  was  considered  by  Pratt  (1970) 
and  results  primarily  from  the  inability  of  the  inlet  jet  to  follow 
the  initial  flame  tube  wall  contour.  This  behavior  leaves  a  "dead" 
zone  or  recirculation  region  adjacent  to  the  flame  tube  wall  from 
which  fluid  is  entrained  by  the  inlet  jet,  thus  establishing  a 
recirculating  flow.  His  predictions  for  this  situation  were  discussed 
previously.  This  mode  of  reverse  flow  would  be  expected  to  diminish 
with  increasing  swirl  of  the  inlet  jet  since  the  jet  will  thus  spread 
faster  and  reach  the  flame  tube  wall  sooner  reducing  the  available 
recirculation  volume.  Decreasing  the  abruptness  of  the  initial  flame 
tube  wall  contour  and  increasing  the  ratio  of  swirler  diameter  to 
Initial  flame  tube  diameter  will  also  diminish  the  magnitude  of  this 
recirculation  mode.  Therefore,  it  is  reasonable  ':o  neglect  this  mode 
except  for  combustors  with  low  swirl  (S  <  1),  abrupt  initial  flame 
tube  contours,  and  small  inlet  jets. 

The  second  mode  of  recirculation  which  occurs  at  high  swirl  can 
assume  two  forms:  the  first  a  reverse  flow  zone  along  the  centerline 
of  the  combustor  and  the  second  an  annular  recirculation  zone  with 
flow  occurring  in  the  normal  direction  in  the  center.  The  first  form 
of  recirculation  was  observed  experimentally  by  Gradon  and  Hiller  (1967), 

Cooper  and  Marek  (1965),  and  Chigier  and  Chervinsky  (1967).  Chigier 
and  Chervinsky  covered  the  low  swirl  range  up  to  S  =  .64  which  marked 
the  onset  of  centerline  reverse  flow;  they  were  able  to  correlate  their 
axial  velocity  profiles  as  error  curves.  Pratt  (197C)  outlined  an 
analytical  technique  for  determining  the  tracer  response  characteristics 
for  this  case,  but  no  empirical  or  analytical  solution  has  been  found. 

The  central  core  recirculation  is  usually  evidenced  in  gas  turbine 
combustors. 

The  annular  type  of  recirculation  was  observed  experimentally  by 
Hiett  and  Powell  (1962).  This  type  of  recirculation  (Fig.  23)  appears 
to  result  from  the  influence  of  the  fuel  injector  flow  which  penetrates 
the  center  of  the  reverse  flow  region  establishing  normal  flow  along 
the  centerline.  When  air-atomizing  fuel  injectors  are  used  with  high 
air  flows  this  type  of  flow  is  likely. 

The  third  mode  of  reverse  flow  is  distinct  from  the  other  two  which 


recirculate  hot,  reacting  gases.  In  this  mode  fresh  secondary  combus¬ 
tion  air  is  forced  upstream  as  a  result  of  penetration  jet  interaction. 
The  fraction  of  secondary  air  flow  directed  into  the  primary  zone  was 
given  empirically  by  Rosenthal  (Anon.  1968)  as: 


mS,P  _ 
mS 


0.5  sin  $  ( 


Ts  1/2 


(4-1) 


where  <j>  is  the  penetration  jet  angle,  is  the  jet  temperature,  and 
Tp  is  the  mean  primary  zone  temperature. 

Verduzio  and  Campanaro  (1969)  performed  a  theoretical  analysis  of 
the  secondary  jet  recirculation  in  tubular  combustors  by  solving  the 
axial  momentum  equation  for  jet  impingement.  Their  equations  are  too 
numerous  to  be  presented  here,  but  the  analysis  showed  that  the  fraction 
of  recirculated  flow  is  dependent  on  the  ratios  of  jet  and  primary  mass 
flows  and  densities,  the  primary  zone  equivalence  ratio,  and  the  ratio 
of  secondary  hole  diameter  to  flame  tube  diameter.  They  were  also  able 
to  show  good  agreement  with  experimental  data.  The  flow  pattern  of 
Fig.  23  shows  this  type  of  reverse  flow,  and  in  general  one  would  expect 
it  to  be  significant  in  gas  turbine  combustors. 

The  primary  zone  of  gas  turbine  combustors  will  be  dominated  by 
the  latter  two  recirculation  modes.  Of  these  two  modes  the  secondary 
jet  impingement  reverse  flow  is  fairly  well  understood  both  analytically 
and  empirically,  while  the  central  recirculation  for  realistic  values 
of  swirl  is  poorly  understood.  The  lack  of  good  definition  of  these 
primary  zone  flows  will  necessarily  hamper  the  primary  zone  modelling 
effort. 

The  secondary  combustion  zones  do  not  exhibit  any  reverse  flow 
regions  (or  other  peculiarities)  and  the  problem  here  reduces  to  model¬ 
ling  the  turbulent  mixing  of  the  secondary  and  dilution  air  jets  and 
cooling  jets  with  the  reacting  main  flow  inside  the  flame  tube.  Tur¬ 
bulent  jet  mixing  lias  been  studied  predominantly  by  empirical  techniques. 
A  good  summary  of  the  available  empirical  results  is  presented  in  the 
work  of  N.R.E.C.  for  NASA  Lewis  Research  Center  (Anon.,  1968),  and  we 
will  utilize  their  results  for  computing  the  rate  of  air  addition  in 
the  secondary  zones.  A  typical  analytical  jet  mixing  approach  is 
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given  by  Hubble  (1967). 

We  do  have  available  a  conputer  program  from  NASA  Lewis  Research 
Center  (Anon.  1968)  capable  of  accurately  predicting  the  air  flow  and 
pressure  distributions  for  a  gas  turbine  combustor  (Tacina  and  Grobman, 
1969).  This  program  has  three  major  limitations  which  restrict  its 
utility  from  a  practical  viewpoint: 

1.  The  primary  zone  is  treated  as  a  single  stirred  reactor; 

2.  The  flow  Inside  the  flame  tube  is  considered  as  being  one-dimen¬ 
sional  ; 

3.  The  fraction  of  fuel  consumed  as  a  function  of  axial  position  in 
the  combustor  must  be  input. 

If  our  combustion  model  can  predict  the  fuel  consumption  function 
accurately,  the  error  incurred  in  the  airflow  and  pressure  distributions 
from  the  first  assunption  will  be  reduced.  It  is  possible  that  an 
inter-program  iteration  will  be  required  if  the  predicted  airflow  dis¬ 
tribution  is  heavily  dependent  on  the  fuel  consumption  function,  but 
presently  this  is  not  thought  to  be  necessary.  Since  this  program  can 
compute  the  air  addition  by  various  jet  mixing  models,  a  realistic 
estimate  of  the  air  available  for  combustion  can  be  incorporated  into 
the  combustion  model.  A  major  difficulty  in  this  program  is  that  it 
considers  only  the  third  recirculation  mode  and  completely  neglects 
the  primary  zone  flow  structure. 

We  have  now  considered  tne  major  flow  characteristics  found  in 
gas  turbine  combustors,  at  least  qualitatively.  The  accurate  airflow 
and  pressure  distributions  predicted  by  the  NASA  combustor  analysis 
program  can  be  used  as  input  to  the  combustion  model.  It  remains  to 
formulate  an  overall  combustion  model  incorporating  a  realistic  primary 
zone  model  and  integrating  the  predicted  airflow  distribution  with  the 
reaction  kinetics  in  the  secondary  combustion  zones. 

B.  A  refined  analytical  gas  turbine  combustor  model 

It  has  been  previously  postulated  that  as  a  direct  result  of  the 
high  turbulence  level  present  in  gas  turbine  combustors  and  the  cor¬ 
respondingly  high  mixing  rate  that  certain  zones  inside  the  flame  tube 
could  approach  the  stirred  reactor  in  behavior.  By  considering  the 
overall  features  of  the  flow  pattern  these  zones  were  connected  to 


form  an  overall  combustor  node!  (cf.  Hammond  and  Mellor,  1970a, b).  The 
basic  overall  model  has  subsequently  bee  i  refined  to  reflect  the  results 
obtained  In  the  previous  chapters,  and  this  refined  model  will  be  pre¬ 
sented  here.  The  refined  model  Is  shown  schematically  In  Fig.  24. 

The  primary  zone  model  remains  much  the  same;  however  an  additional 
fuel  Inlet  Into  the  recirculation  reactor  Is  considered  after  Hottel 
et  al .  (1958).  The  additional  fuel  distribution  point  permits  studying 
the  effects  of  the  initial  reactant  distribution.  Addition  of  fuel  into 
both  zones  more  closely  corresponds  to  the  physical  situation  in  which 
fuel  vapor  Is  entrained  by  both  swirler  air  and  hot  recirculating  gases. 
The  two  reactor  structure  is  retained  for  simplicity  with  the  main 
flow  reactor  accounting  for  reaction  in  flow  of  the  normal  direction, 
and  the  single  recirculation  reactor  must  account  for  the  reaction  in 
all  three  possible  recirculation  flow  zones.  Note  that  although  this 
model  can  predict  a  wide  variation  in  primary  zone  residence  times  and 
some  variation  in  initial  reactant  distributions,  it  lacks  the  ability 
to  predict  a  characteristic  delay  time  as  observed  by  Beer  and  lee  (1965). 
This  deficiency  can  be  easily  rectified  by  the  Inclusion  of  a  series 
plug  flow  reactor  in  the  primary  zone  or  by  the  recirculating  plug  flow 
pair  of  Morgan  (1967)  shown  in  Fig.  25.  The  ability  of  the  recircul¬ 
ation  pair  model  to  predict  flame  stabilization  phenomena  has  been 
discussed  previously  (Hammond  and  Mel  lor,  1970a, b;  Hottel  et  al.,  1958). 

Our  model  of  the  secondary  combustion  zones  has  also  been  slightly 
altered.  Each  of  the  “plug  flow  reactors"  is  now  modelled  as  a  sequence 
of  perfectly  stirred  reactors  of  finite  and  equal  volume.  Altering 
the  reactor  volumes  gives  a  variation  In  the  residence  time  distribution 
as  discussed  in  the  tanks-in-series  model,  since  the  standard 
deviation  o f  the  residence  time  distribution  depends  on  the  number  of 
reactors  incorporated  in  a  given  volume  (cf.  Eqn.  3-7).  Each  individual 
perfectly  stirred  reactor  now  has  its  own  discrete  amount  of  air  and 
fuel  addition;  therefore  the  predicted  mixing  behavior  can  be  altered 
by  changing  the  fuel  distribution  function.  The  semi -continuous  air 
additions  permit  closer  approximations  of  the  airflow  distribution  com¬ 
puted  using  turbulent  jet  mixing  phenomena.  This  is  particulari ly  im¬ 
portant  when  considering  film-cooled  combustors.  By  varying  the  volumes 
of  the  perfectly  stirred  reactors,  the  spread  (or  deviation)  of  the 


Fig.  25  Recirculation  model  using  a  plug  flow 
reactor  pair  (after  Morgan,  1967). 
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residence  time  distribution  can  be  changed  as  a  function  of  axial 
position. 

The  current  overall  mc^el  has  several  Intrinsic  merits.  Firstly, 

It  Is  simple  enough  to  facilitate  computer  solution  without  consuming 
vast  amounts  of  time.  Secondly,  both  mixing  and  kinetics  are  considered. 
Thirdly,  an  exact  correspondence  exists  between  the  flew  regions  in 
the  model  and  those  in  the  physical  combustor;  therefore,  the  internal 
composition  and  temperature  patterns  will  be  approximately  predicted. 

This  ability  to  approximate  the  internal  structure  of  the  reacting 
flow  separates  the  model  from  those  relying  heavily  on  exit  tracer  re¬ 
sponse  functions  to  predict  overall  performance  and  which  neglect  in¬ 
ternal  structure  entirely.  In  fact,  to  obtain  more  exact  and  detailed 
results  the  coupled  turbulent  governing  equations  must  be  integrated. 

Having  presented  what  feel  to  be  a  realistic  overall  combustor 
model  let  us  examima  the  c  rrent  nodels  of  other  investigators  before 
presenting  the  predictions  of  our  own  model. 

C.  Other  recently  proposed  combustor  models 

There  exist  only  two  other  complete  analytical  models  for  gas  tur¬ 
bine  combustors.  The  first  was  proposed  by  Heywood  et  al .  (1970)  for 
use  in  predicting  air  pollutant  emissions.  They  modelled  the  primary 
combustion  zone  as  a  single  plug  flow  reactor  and  the  secondary  com¬ 
bustion  zones  as  a  one-dimensional  equilibrium  flow  with  continuous  air 
addition.  Although  this  model  gave  fair  agreement  for  nitric  oxide 
emissions,  because  it  assumes  perfect  mixing  in  the  radial  direction 
In  the  primary  zone  and  thus  neglects  the  recirculation  flow,  it  cannot 
predict  combustor  blowout  and  flame  stability.  The  use  of  C-H-0  equITI- 
brlum  neglects  the  kinetic  nature  of  nitric  oxide  formation  and  can 
lead  to  serious  error  as  shown  by  Bowman  (1970). 

The  other  complete  gas  turbine  combustor  model  was  developed  by 
Fletcher  and  Heywood  (1971)  to  predict  the  nitric  oxide  emission  levels 
of  aircraft  gas  turbine  combustors.  The  primary  zone  is  modelled  as 
a  single  segregated  perfectly  stirred  reactor  with  an  assumed  initial 
normal  distribution  of  mixture  ratio.  The  effects  of  such  assumptions 
have  been  discussed  in  the  section  on  segregation.  An  arbitrary 
fraction  of  the  total  fuel  flow  is  assumed  to  burn  in  the  primary  zone. 
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the  rest  being  carried  downstream.  The  primary  zone  Is  followed  by  an 
intermediate  zone  which  Is  modelled  by  forcing  the  standard  deviation 
of  the  mixture  ratio  distribution  to  zero  over  an  arbitrary  zone  length. 
This  corresponds  to  the  intermixing  of  eddies  in  a  segregated  plug  flow, 
and  at  the  exit  of  the  zone  the  flow  is  micromixed  plug.  The  final  or 
dilution  zone  is  then  modelled  as  micromixed  plug  flow  with  air  addition. 

One  major  limitation  of  this  model  is  that  it  holds  the  carbon- 
hydrogen-oxygen  chemistry  in  equilibrium  and  treats  only  the  nitric 
oxide  formatiof.  kinetically,  since  it  is  assumed  that  the  entire  pro¬ 
cess  is  mixing  limited.  As  mentioned  previously,  the  work  of  Bowman 
(1970)  discredits  the  equilibrium  assumption.  Also,  this  model  contains 
a  nunber  of  arbitrary  parameters  which  must  be  assigned  values  intuitively 
(as  does  our  model  to  a  lesser  degree). 

The  predicted  nitric  oxide  emission  levels  exhibited  the  correct 
trends  and  reasonable  values.  The  strong  point  of  this  model  (and  our 
refined  model)  is  its  ability  to  examine  the  effect  of  segregation  of 
primary  zone  mixture  ratio.  The  initial  calculations  used  only  the 
primary  and  dilution  zones,  but  showed  an  interesting  variation  in 
nitric  oxide  emissions  with  initial  reactant  distribution  (Fig.  26). 

This  result  indicates  that  the  primary  zone  mixing  should  receive 
attention  in  our  future  work. 

We  now  see  that  most  investigators  agree  on  the  basic  premise  of 
using  combinations  of  perfectly  stirred  and  plug  flow  reactors  to  model 
gas  turbine  combustors.  However,  there  exist  differences  in  opinion 
as  to  the  correct  techniques  for  handling  segregation,  mixing,  and 
kinetics.  The  true  test  in  each  case  will  be  the  experimental  confir¬ 
mation  of  the  analytical  predictions.  For  this  reason  we  now  turn  to 
a  discussion  of  our  predicted  results. 
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CHAPTER  V 

PREDICTIONS  OF  THE  ANALYTICAL  GAS 
TURBINE  COMBUSTOR  MODEL 


As  a  result  of  difficulties  with  the  NASA  combustor  program 
(Anon.,  1968)  only  their  test  case  air  flow  and  pressure  distributions 
have  been  available.  He  have,  however,  analyzed  the  combustion  pro¬ 
cess  in  this  test  case  combustor  (denoted  No.  3)  using  the  overall 
combustor  analysis  program  (unfortunately  the  NASA  Test  Case  Combustor 
No.  3  Is  not  known  to  correspond  to  ary  actual  combustor).  Although 
this  Is  an  annular  configuration  our  analysis  applies  equally  well. 

To  reduce  computation  time,  the  combustor  was  modelled  using  single 
perfectly  stirred  reactors  for  all  secondary  combustion  zones.  The 
heavy  hydrocarbon  fuel  originally  used  was  replaced  with  an  equivalent 
flow  of  propane,  all  of  which  was  Injected  into  Vj.  Arbitrarily,  a 
recirculation  ratio  R  ■  0.5  and  a  primary  zone  volune  split  =  1.0 

were  assured.  The  configuration,  airflow,  ard  pressure  distributions 
taken  from  the  NASA  output  and  used  as  Input  to  the  modelling  program 
are  shown  in  Figs.  27-29.  The  other  operating  conditions  are: 


air  mass  flow  rate 

Inlet  air  temperature 

fuel  mass  flow  rate 

Inlet  fuel  temperature 

mean  primary  zone  equivalence  ratio 


■  4358  g/sec 
-  895  °K 

*  40  g/sec 

*  293  °K 
»  1.34 


The  details  of  the  gas  turbine  combustor  analysis  program  are  given  in 
Appendix  C.  We  will  now  discuss  the  predictions  of  the  program  both 
from  a  performance  and  from  a  pollution  viewpoint. 

A.  Performs  nee  predictions 

The  temperature,  completeness  of  combustion,  and  the  mole  fractions 


50.00 


0-001 


of  all  species  have  been  plotted  versus  axial  position  In  the  combustor 
these  results  are  shown  In  Figs.  30-42.  Those  most  closely  related  to 
combustion  performance  are  the  temperature,  completeness  of  combustion, 
and  carbon  monoxide  mole  fraction  profiles  (Figs.  30,31,  and  41). 

The  temperature  profile  peaks  In  the  primary  zone  and  decreases 
downstream  as  expected.  The  slightly  reduced  teirperature  in  the  recir¬ 
culation  zone  results  from  the  Increased  dissociation  in  this  region 
Indicating  that  the  C-H-0  reactions  in  the  recirculation  zone  have  had 
time  to  more  closely  approach  equilibrium  than  those  in  the  primary 
zone.  This  dissociation  (particularlly  of  carbon  dioxide)  also  reduces 
the  completeness  of  combustion  in  the  recirculation  zone.  Carbon  mon¬ 
oxide  is  shown  to  peak  in  the  recirculation  zone  since  the  high  primary 
zone  equivalence  ratio  favors  the  formation  of  that  species  over  carbon 
dioxide.  The  primary  zone  also  contains  the  maxima  of  all  species 
having  a  highly  reactive  nature  (free  radicals),  most  notably  the  oxygen 
atom  and  the  hydroxyl  radical. 

Considering  results  of  an  overall  nature,  we  see  that  combustion 
is  essentially  complete  in  roughly  half  the  length  of  the  combustor. 
This  indicates  that  a  considerable  length  reduction  could  be  realized 
for  this  combustor  without  significant  performance  degradation  if  a 
turbine  inlet  temperature  of  1500°K  could  be  tolerated.  The  final  com¬ 
bustion  efficiency  is  essentially  1002  in  correspondence  with  values 
observed  in  well-designed  systems.  Note  that  the  predicted  exit  temper¬ 
ature  of  about  1200°K  is  physically  reasonable. 

B.  Pollutant  emissions 

The  primary  pollutant  species  of  interest  are  carbon  monoxide  and 
nitric  oxide  (hydrocarbon  emissions  are  presently  not  predicted  because 
of  the  assumed  Infinitely  fast  and  complete  partial  oxidation  to  CO). 

The  carbon  monoxide  profile  indicates  that  a  negligible  amount  of  this 
species  will  be  emitted.  The  exit  carbon  monoxide  concentration  of  19 
ppm  is  in  agreement  with  typical  values  for  practical  systems  (Cornelius 
and  Wade,  1970). 

Examining  the  predicted  nitric  oxide  concentration  profile  reveals 
two  interesting  points.  Firstly,  the  nitric  oxide  level  peaks  in  the 
recirculation  zone  in  complete  agreement  with  the  concept  that  the 
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FIGURE  31.  COMPLETENESS  OF  COMBUSTION  PREDICTED  FO! 
TEST  CASE  COMBUSTOR  NO.  3  (RECIRC.  RRTI0=0  .SWOLUME 
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amount  of  nitric  oxide  formed  depends  on  the  length  of  time  spent  at 

peak  temperatures.  Secondly,  the  exit  nitric  oxide  concentration  of 

142  ppm  Is  within  the  range  of  values  observed  In  combustors  (Cornelius  i 

and  Wade,  1970).  ] 

k 

We  can  make  no  experimental  comparisons  for  the  NASA  Test  Com-  j 

bustor  No.  3  other  than  the  generalizations  mentioned  above  since  this  j 

combustor  Is  not  a  physical  unit.  Indeed  for  this  reason  no  time-con-  j 

s unlng  parametric  studies  were  performed.  j 

■  .  1 

I 


87 


CHAPTER  VI 

CONCLUSIONS  AND  FUTURE  EFFORTS 

A  complete  model  for  the  analysis  of  all  types  of  gas  turbine 
combustors  has  been  developed.  It  utilizes  the  fluid  mechanical  portion 
of  the  experimentally  verified  NASA  combustor  analysis  program  (Anon., 
1968)  to  compute  required  distributions  of  airflow  and  static  pressure. 
The  combustion  model  itself  consists  of  an  assemblage  of  perfectly 
stirred  reactors.  Inside  each  perfectly  stirred  reactor  a  fully  kinetic 
combustion  mechanism  for  a  general  hydrocarbon  is  used  to  compute  per¬ 
formance.  The  level  of  nitric  oxide  formation  is  also  kinetically  com¬ 
puted.  Upon  consideration  of  the  physical  flows  inside  the  combustor, 
the  vol unes  of  these  perfectly  stirred  reactors,  the  flows  intercon¬ 
necting  them,  and  the  fuel  and  air  addition  points  were  arranged  to 
allow  a  reasonable  approximation  of  the  overall  mixing  behavior  in 
terms  of  a  minimum  number  of  physically  significant  parameters.  The 
predictions  of  the  model  yield  approximate  internal  profiles  of  tem¬ 
perature  and  species  concentrations  which  Indicate  performance  and  can 
be  directly  employed  in  a  heat  transfer  analysis  to  predict  flame  tube 
wall  temperatures  (Owens  and  Mel  lor,  1971). 

The  kinetic  analysis  for  general  hydrocarbon  combustion  is  based 
on  the  postulate  that  the  combustion  process  of  any  hydrocarbon  can  be 
separated  into  three  component  processes: 

1.  a  rapid  partial  oxidation  of  the  parent  hydrocarbon  to  carbon 
monoxide  and  water; 

2.  a  single  rate  limiting  step  of  carbon  monoxide  destruction  (Reaction 

CM]); 

3.  a  sequence  of  free  radical  production  reactions  not  involving  carbon; 

where  step  1  is  assumed  to  have  infinite  rate  compared  to  step  2.  It 
was  deduced  that  the  accuracy  of  these  simplifications  should  improve 
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as  the  complexity  of  the  hydrocarbon  increases. 

The  Zeldovich  mechanism  was  taken  to  be  responsible  for  the  majority 
of  nitric  oxide  formation,  but  the  three-body  association  reaction  of 
monatomic  nitrogen  and  oxygen  was  included  since  it  is  known  to  be  of 
increasing  significance  at  high  temperatures. 

A  study  of  conventional  chemical  reactor  theory  revealed  techniques 
for  modelling  mixing  by  proper  combinations  of  perfectly  stirred  reactors. 
It  was  found  that  an  arbitrary  residence  time  distribution  can  be 
approximated  by  combinations  of  plug  flow  reactors  and  perfectly  stirred 
reactors  with  various  interconnecting  flows.  Investigation  of  segrega¬ 
tion  showed  that  the  residence  time  distribution  alone  is  not  sufficient 
to  predict  behavior  due  to  the  eddy  dominated  nature  of  the  turbulent 
mixing  processes.  In  particular,  if  the  fluid  segregates  the  initial 
distribution  of  reactants  over  the  eddies,  the  rate  of  eddy  intermixing, 
and  the  eddy  residence  time  distribution  must  be  determined.  By  con¬ 
sidering  various  approximations  of  these  functions  a  number  of  stirred 
reactors  were  interconnected  to  approximate  the  mixing  processes. 

Tha  resulting  overall  combustor  model  was  seen  to  posses*-,  several 
intrinsic  merits,  namely: 

1.  it  is  simple  enough  to  facilitate  solutions; 

2.  both  mixing  and  kinetics  are  considered; 

3.  there  is  an  exact  correspondence  between  the  flow  zones  in  the 
model  and  in  a  combustor. 

Therefore,  it  is  inherently  superior  to  those  models  which  neglect 
mixing  and/or  kinetics,  or  those  which  rely  heavily  on  exit  tracer  re¬ 
sponse  functions  and  neglect  internal  structure. 

It  is  presumptuous  to  make  conclusions  concerning  the  predictions 
of  this  model  in  view  of  the  lack  of  experimental  data.  However,  trend 
predictions  inside  the  combustor  are  seen  to  be  accurate,  and  the  values 
of  the  variables  at  combustor  exit  fall  in  the  range  of  values  observed 
In  experimental  systems. 

The  major  region  of  future  refinements  will  be  the  primary  zone 
model.  A  need  exists  for  good  work  on  the  three  recirculation  modes 
and  the  exact  influence  of  each  on  mixing  and  reaction.  Relatively 
little  is  known  concerning  the  internal  structure  of  the  initial  vapori- 
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zation  zone.  In  particular  the  droplet  distribution  and  accompanying 
evaporation  rates.  The  question  concerning  the  existence  of  hetero¬ 
geneous  combustion  is  yet  unresolved.  The  modes  of  entrainment  of 
fuel  vapor  from  the  zone  will  define  the  initial  mixture  ratio  segre¬ 
gation  In  the  primary  combustion  zone,  but  this  is  poorly  understood. 

Future  work  using  the  current  analytical  model  will  include  ex¬ 
tensive  parametric  studies  of  the  Allison  J-33  combustor  and  the  Boeing 
502E  combustor  for  several  operating  points.  These  studies  should 
elucidate  the  erfects  of  the  various  mixing  parameters.  These  predic¬ 
tions  will  then  be  compared  with  experimental  data  on  these  combustors 
for  the  same  operating  points.  The  experimental  facility  and  technique 
is  described  fully  by  Anderson  and  Mellor  (1971).  Comparison  will  be 
made  on  the  basis  of  the  internal  profiles  of  temperature,  carbon 
monoxide  concentration  and  nitric  oxide  concentration.  It  is  hoped 
that  further  refinements  of  the  combustor  modelling  program  will  become 
apparent  as  a  result  of  these  comparisons. 
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APPENDIX  A 

DESCRIPTION  OF  THE  KINETIC  PERFECTLY  STIRRED 
REACTOR  PROGRAM  FOR  HYDROCARBON  COMBUSTION 

The  function  of  this  computer  program  is  to  predict  the  perform¬ 
ance  of  a  single  stirred  reactor  by  solving  the  coupled  set  of  equa¬ 
tions  (2-5),  (2-6),  and  (2-7).  The  current  program  is  an  extensive 
revision  of  the  program  developed  by  Jones  and  Prothero  (1968)  which 
solved  the  set  of  equations  (2-5)  and  (2-6),  i.e.  solved  for  the 
concentrations  for  a  specified  reactor  temperature.  This  formulation 
was  unacceptable  for  use  in  the  overall  model;  thus  we  included  Eqn. 
(2-7)  in  the  set  allowing  solutions  for  a  specified  inlet  temper¬ 
ature.  Curve  fits  for  the  thermochemical  properties  were  abstracted 
from  the  work  of  Prothero  (1969).  The  original  kinetics  data  were 
altered  in  accordance  with  Table  9.  A  brief  description  of  the  main 
program  and  each  subprogram  follows  (Jones,  1969). 

THORTNA 

THORTNA  is  the  main  control  program  which  functions  mainly  as  an 
input/output  routine.  Optionally,  two  levels  of  printed  output  and 
also  punched  card  output  will  be  produced. 

INPTHD 

INPTHD  is  a  library  data  input  routine.  It  reads  a  standard 
thermochemical  data  package  from  magnetic  tape  at  the  start  of  each 
program  run.  This  data  package  contains  values  of  the  specific  heat 
curve  fit  coefficients,  the  molecular  weights ,  enthalpies  and  entropies 
of  formation,  and  several  arrays  which  define  the  various  reactions. 
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SLCTNfl 

SLCTNA  sets  up  the  mechanism  and  rate  constant  matrices  for  the 
reactions  chosen  by  selecting  data  from  the  library  data  package. 

TSOLVE 

TSOLVE  performs  the  solution  of  the  energy  equation  by  matching 
inlet  and  exit  specific  enthalpies.  It  uses  a  golden  section  tech¬ 
nique  and  calls  GASMIX  once  per  iteration. 

GASHI X 

Given  an  estimated  reactor  temperature  from  TSOLVE,  GASMIX  solves 
for  the  density  and  species  concentrations  existing  in  the  reactor.  A 
multi-dimensional  extension  of  Newton's  method  is  used.  Both  FUNCTN 
and  DERIVS  are  called. 

FUNCTN 

FUNCTN  evaluates  the  residuals  of  the  species  conservation  equa¬ 
tions  which  are  written  in  the  form  F  (a^ ,  p)  =  0. 

DERIVS 

DLRIVS  evaluates  the  Jacobian  matrix  of  the  partial  derivatives 
of  the  species  conservation  equations. 

SRENTHA  -  SRH298A 

These  function  subprograms  operate  together  to  evaluate  the  molar 
enthalpies  of  a  particular  species  as  a  function  of  temperature. 

HATINV  -  MXPERM 

These  two  subroutines  are  called  by  GASMIX  to  invert  the  Jacobian 
matrix  of  partial  derivatives. 

KDATA 

KDATA  is  a  block  data  subprogram  which  stores  the  values  of  the 
parameters  which  define  the  various  rate  constants. 

To  use  this  program  for  general  hydrocarbon  combustion  the  parent 
fuel  must  be  converted  to  its  equivalent  in  carbon  monoxide  and  water. 
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and  an  appropriate  heat  transfer  must  be  computed.  This  process  is 
accomplished  by  the  following  global  reaction: 

h  f 

4>C^Hj^  +  (a  +  ij-)  O2  .  $aC0  +  <i> ( b/2 )  H20  +  £a(l-<j>/2)  + 

(l-<)>)b/4]  02  +  Q 

where  the  right  hand  side  gives  the  true  inlet  mixture  to  the  com¬ 
bustor,  and  Q  is  the  heat  transfer  to  the  reactor  when  expressed  as 
cal/gm  of  flow.  Here  $is  the  equivalence  ratio  of  the  original  hydro 
carbon  mixture. 

This  program  is  completely  operational  and  detailed  information 
will  be  supplied  upon  request.  Typical  computation  times  on  a 
COC  6500  machine  range  from  30  to  150  sec  per  case  depending  on  the 
severity  of  the  conditions  imposed. 
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APPENDIX  B 

DESCRIPTION  OF  THE  EQUILIBRIUM  REACTOR  PROGRAM 

This  program  was  designed  to  give  an  upper  limit  on  the  pre¬ 
dictions  of  the  kinetic  program  (Appendix  A).  As  a  result  it  utilizes 
the  same  library  data  package  but  processes  the  information  in  a 
slightly  different  manner.  The  program  is  an  inhouse  product  and  has 
not  been  previously  presented.  The  numerical  technique  is  rather 
unique:  starting  from  some  initial  estimate  of  the  composition,  the 
current  composition  is  used  to  evaluate  all  of  the  equilibrium  con¬ 
stants  and  a  temperature  from  the  energy  equation.  A  temperature  is 
then  calculated  from  each  of  the  equilibrium  constants.  The  devi¬ 
ation  of  each  of  these  temperatures  from  the  average  is  used  to 
increment  the  composition,  and  the  process  repeats  until  some  con¬ 
vergence  criterion  is  met.  Solutions  can  be  performed  for  fixed 
reactor  or  inlet  temperature.  Again  each  subprogram  will  be  briefly 
described. 

EQRTR 

EQRTR  is  the  main  program  which  serves  only  as  an  input/output 
routine.  Optionally  both  printed  and  punched  card  output  can  be  pro¬ 
duced.  *, 

GSOLV 

GSOLV  is  an  organizational  routine  which  controls  the  initial¬ 
ization  and  calling  sequence  of  the  operational  routines. 

EQDTIN 

EQDTIN  is  the  library  data  input  routine  for  this  program. 

NEWNI 

NEWNI  computes  the  composition  increments  mentioned  in  discussing 
the  numerical  technique. 
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PRISOL 

PRISOL  is  used  by  NEWNI  in  computing  the  composition  increments. 
It  solves  the  set  of  atomic  species  conservation  equations. 

TEMPKP 

TEM’KP  evaluates  a  temperature  from  each  of  the  equilibrium 
constant  values. 

TENERG 

TENERG  evaluates  a  temperature  by  solving  the  energy  equation. 

Typical  computation  times  for  this  program  on  the  CDC  6500 
machine  range  from  20  to  60  sec  depending  heavily  on  the  accuracy 
of  the  optionally  supplied  initial  estimate.  We  will  supply,  on 
request,  a  more  complete  description  of  this  program. 
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APPENDIX  C 

DESCRIPTION  OK  THE  GAS  TURBINE  COf-BUSTOR  ANALYSIS  PROGRAM 

This  program  performs  the  analysis  of  a  complete  gas  turbine 
combustor  based  on  the  overall  model.  The  basic  comoutational  device 
is  the  kinetic  perfectly  stirred  reactor  program  described  in  Appendix 
A.  Certain  inputs  to  the  program  are  required  and  they  will  be 
summarized  here: 

1.  the  combustor  geometry; 

2.  the  predicted  airflow  and  pressure  distributions  from  the 
NASA  program; 

3.  an  estimated  fuel  distribution; 

4.  the  total  fuel  flow  and  temperature; 

5.  the  total  airflow  and  temperature; 

6.  estimates  of  the  primary  zone  recirculation  ratio  and 
volume  ratio; 

7.  the  species  and  reactions  to  be  considered. 

The  same  thermochemical  data  package  is  utilized.  To  compact  the 
program  as  much  as  possible  it  was  written  in  overlayed  form  (cf.  Fig. 
43).  To  conserve  computer  time  an  internal  abort/restart  sequence  was 
written  into  the  program.  The  abort  sequence  functions  if  the  program 
detects  an  error  (insufficient  time,  incorrect  input,  non-convergence, 
etc.)  and  at  that  point  all  internal  information  is  dumped  to  magnetic 
tape.  Using  the  restart  sequence  the  program  can  be  started  at  any 
point  either  for  error  correction  or  parametric  studies.  For  instance, 
when  running  two  cases  which  differ  only  in  secondary  air  addition, 
the  program  could  be  started  at  the  point  where  the  distribution  first 
diff.rs,  thereby  using  a  portion  of  the  previous  calculations.  This 
feature  will  be  useful  when  changing  the  secondary  combustion  zones 
while  the  primary  zone  description  remains  fixed.  A  description  of 
each  routine's  function  follows. 
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GTGO'-Sl 

This  is  the  resident  main  control  program;  it  directs  the  calling 
sequence  of  the  other  overlays  and  computes  the  total  time  consumed 
by  each.  It  also  contains  all  the  input/output  buffers. 

ABORT 

As  the  name  implies  this  is  the  abort  sequence  routine.  If  an 
error  is  detected,  this  routine  is  called,  and  it  dumps  all  internal 
information  to  magnetic  tape.  An  appropriate  error  message  is  issued 
to  the  output  device.  ABORT  is  also  called  at  the  completion  of  com¬ 
putations  to  save  the  internal  information  for  future  runs. 

INPUT2 

INPUT2  is  the  main  program  of  the  library  data  input  overlay.  It 
performs  two  calls,  one  on  INPTHD  and  one  on  TIMEL.  The  function  of 
INPTHD  was  described  in  Appendix  A. 


TIMEL 


TIMEL  is  an  assembly  language  coded  routine  v'  evaluates  the 
total  time  available  to  the  program  f<*r  -.caution.  This  time  is  used 
in  the  abort  sequence.  .... 

inputi 


ivr.itl  is  the  main  program  of  the  case  data  input  overlay.  It 
^-/.V'CCds  the  previously  mentioned  input  data,  and  sets  up  the  combustor 
model  by  assigning  volumes,  flows,  and  pressures  to  all  of  the  reactors 


QTFG 

QTFG  is  used  by  INPUTI  to  integrate  the  combustor  dimensions  in 
order  to  compute  the  volume. 

ATSM  -  ACFI 

These  two  subroutines  are  used  by  INPUT1  and  OUTPT  to  interpolate 
values  for  the  various  reactor  parameters  from  the  tabulated  distri¬ 
bution  functions. 


THORTNB 

THORTNB  is  the  main  program  of  the  computational  overlay.  It 
controls  the  processing  in  that  overlay. 

FLOMIX 


FLOMIX  is  used  to  initialize  the  input  conditions  for  each  re¬ 
actor.  It  performs  the  partial  oxidation  calculation  discussed  in 
Appendix  A. 


>t;ci;iRC 

RECIRC  solves  for  the  performance  of  the  primary  zone  model. 


REACTR 


REACTR  performs  the  actual  solution  for  each  reactor.  It  is  a 
slight  revision  of  TSOLVE  described  in  Appendix  A. 

TRPRNT 

TRPRNT  serves  as  an  intermediate  output  routine  for  debugging 
purposes. 

The  remainder  of  the  computational  routines  were  described  in 
Appendix  A. 


OUTPT 

OUTPT  prints  out  the  solution  predicted  by  the  computational  over 
lay.  In  addition  it  can  plot  the  input  distributions  and  output 
profiles  as  functions  of  axial  position  and/or  mean  residence  time. 

The  exact  configuration  of  this  program  is  in  a  state  of  flux. 

For  the  test  case  presented  in  Chapter  V  the  total  computation  time 
was  about  1100  sec  on  the  CDC  6500  machine. 


